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Abstract 

These notes are intended as an introductory course for experimental particle 
physicists interested in the recent developments in astrophysics and cosmo- 
logy. I will describe the standard Big Bang theory of the evolution of the 
universe, with its successes and shortcomings, which will lead to inflation- 
ary cosmology as the paradigm for the origin of the global structure of the 
universe as well as the origin of the spectrum of density perturbations respon- 
sible for structure in our local patch. I will present a review of the very rich 
phenomenology that we have in cosmology today, as well as evidence for the 
observational revolution that this field is going through, which will provide us, 
in the next few years, with an accurate determination of the parameters of our 
standard cosmological model. 

1 GENERAL INTRODUCTION 

Cosmology (from the Greek: kosmos, universe, world, order, and logos, word, theory) is probably the 
most ancient body of knowledge, dating from as far back as the predictions of seasons by early civiliza- 
tions. Yet, until recently, we could only answer to some of its more basic questions with an order of 
magnitude estimate. This poor state of affairs has dramatically changed in the last few years, thanks to 
(what else?) raw data, coming from precise measurements of a wide range of cosmological parameters. 
Furthermore, we are entering a precision era in cosmology, and soon most of our observables will be 
measured with a few percent accuracy. We are truly living in the Golden Age of Cosmology. It is a very 
exciting time and I will try to communicate this enthusiasm to you. 

Important results are coming out almost every month from a large set of experiments, which pro- 
vide crucial information about the universe origin and evolution; so rapidly that these notes will probably 
be outdated before they are in print as a CERN report. In fact, some of the results I mentioned during 
the Summer School have already been improved, specially in the area of the microwave background 
anisotropics. Nevertheless, most of the new data can be interpreted within a coherent framework known 
as the standard cosmological model, based on the Big Bang theory of the universe and the inflationary 
paradigm, which is with us for two decades. I will try to make such a theoretical model accesible to 
young experimental particle physicists with little or no previous knowledge about general relativity and 
curved space-time, but with some knowledge of quantum field theory and the standard model of particle 
physics. 



2 INTRODUCTION TO BIG BANG COSMOLOGY 

Our present understanding of the universe is based upon the successful hot Big Bang theory, which 
explains its evolution from the first fraction of a second to our present age, around 13 bilhon years 
later. This theory rests upon four strong pillars, a theoretical framework based on general relativity, 
as put forward by Albert Einstein ^ and Alexander A. Friedmann ^ in the 1920s, and three robust 
observational facts: First, the expansion of the universe, discovered by Edwin P. Hubble ^ in the 1930s, 
as a recession of galaxies at a speed proportional to their distance from us. Second, the relative abundance 
of light elements, explained by George Gamow [Q] in the 1940s, mainly that of helium, deuterium and 
lithium, which were cooked from the nuclear reactions that took place at around a second to a few minutes 



after the Big Bang, when the universe was a few times hotter than the core of the sun. Third, the cosmic 
microwave background (CMB), the afterglow of the Big Bang, discovered in 1965 by Arno A. Penzias 
and Robert W. Wilson as a very isotropic blackbody radiation at a temperature of about 3 degrees 
Kelvin, emitted when the universe was cold enough to form neutral atoms, and photons decoupled from 
matter, approximately 500,000 years after the Big Bang. Today, these observations are confirmed to 
within a few percent accuracy, and have helped establish the hot Big Bang as the preferred model of the 
universe. 



2.1 Friedmann-Robertson-Walker universes 

Where are we in the universe? During our lectures, of course, we were in Casta Papiemicka, in "the heart 
of Europe", on planet Earth, rotating (8 light-minutes away) around the Sun, an ordinary star 8.5 kpc[] 
from the center of our galaxy, the Milky Way, which is part of the local group, within the Virgo cluster 
of galaxies (of size a few Mpc), itself part of a supercluster (of size ~ 100 Mpc), within the visible 
universe (~ few x 1000 Mpc), most probably a tiny homogeneous patch of the infinite global structure 
of space-time, much beyond our observable universe. 

Cosmology studies the universe as we see it. Due to our inherent inability to experiment with it, 
its origin and evolution has always been prone to wild speculation. However, cosmology was born as a 
science with the advent of general relativity and the realization that the geometry of space-time, and thus 
the general attraction of matter, is determined by the energy content of the universe [^, 

G^v = R^lv - ^9iiuR = SirGTi^iy + Ag^,^ . (1) 

These non-linear equations are simply too difficult to solve without some insight coming from the sym- 
metries of the problem at hand: the universe itself. At the time (1917-1922) the known (observed) 
universe extended a few hundreds of parsecs away, to the galaxies in the local group, Andromeda and 
the Large and Small Magellanic Clouds: The universe looked extremely anisotropic. Nevertheless, both 
Einstein and Friedmann speculated that the most "reasonable" symmetry for the universe at large should 
be homogeneity at all points, and thus isotropy. It was not until the detection, a few decades later, of the 
microwave background by Penzias and Wilson that this important assumption was finally put onto firm 
experimental ground. So, what is the most general metric satisfying homogeneity and isotropy at large 
scales? The Friedmann-Robertson-Walker (FRW) metric, written here in terms of the invariant geodesic 
distance ds^ = g^ydx^^dx^ in four dimensions, = 0, 1, 2, 3, see Ref. [^],| 
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characterized by just two quantities, a scale factor a{t), which determines the physical size of the uni- 
verse, and a constant K, which characterizes the spatial curvature of the universe, 
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Spatially open, flat and closed universes have different geometries. Light geodesies on these universes 
behave differently, and thus could in principle be distinguished observationally, as we shall discuss later. 
Apart from the three-dimensional spatial curvature, we can also compute a four-dimensional space-time 
curvature, 

a /a\2 K 
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*One parallax second (1 pc), parsec for short, corresponds to a distance of about 3.26 light-years or 3 x lO'^* cm. 
am using c — 1 everywhere, unless specified. 



Depending on the dynamics (and thus on the matter/energy content) of the universe, we will have dif- 
ferent possible outcomes of its evolution. The universe may expand for ever, recoUapse in the future or 
approach an asymptotic state in between. 



2.1.1 The expansion of the universe 

In 1929, Edwin P. Hubble observed a redshift in the spectra of distant galaxies, which indicated that they 
were receding from us at a velocity proportional to their distance to us This was correctly interpreted 
as mainly due to the expansion of the universe, that is, to the fact that the scale factor today is larger 
than when the photons were emitted by the observed galaxies. For simplicity, consider the metric of a 
spatially flat universe, ds^ = dt^ — a?{t)dx^ (the generalization of the following argument to curved 
space is straightforward). The scale factor a{t) gives physical size to the spatial coordinates x, and the 
expansion is nothing but a change of scale (of spatial units) with time. Except for peculiar velocities, i.e. 
motion due to the local attraction of matter, galaxies do not move in coordinate space, it is the space-time 
fabric which is stretching between galaxies. Due to this continuous stretching, the observed wavelength 
of photons coming from distant objects is greater than when they were emitted by a factor precisely equal 
to the ratio of scale factors, 

Aobs «0 _ , , 

- — = — = 1 + z, (5) 
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where ag is the present value of the scale factor. Since the universe today is larger than in the past, the 
observed wavelengths will be shifted towards the red, or redshifted, by an amount characterized by z, the 
redshift parameter. 

In the context of a FRW metric, the universe expansion is characterized by a quantity known as 
the Hubble rate of expansion, H{t) = d{t)/a{t), whose value today is denoted by Hq. As I shall deduce 
later, it is possible to compute the relation between the physical distance and the present rate of 
expansion, in terms of the redshift parameter,^ 

HodL = z + ^{l-qo)z^ + 0{z^). (6) 

At small distances from us, i.e. at z ^ 1, we can safely keep only the linear term, and thus the recession 
velocity becomes proportional to the distance from us, v = c z = Hodi, the proportionality constant 
being the Hubble rate, Hq. This expression constitutes the so-called Hubble law, and is spectacularly 
confirmed by a huge range of data, up to distances of hundreds of megaparsecs. In fact, only recently 
measurements from very bright and distant supemovae, at z ~ 1, were obtained, and are beginning to 



probe the second-order term, proportional to the deceleration parameter qo, see Eq. (|22|). I will come 
back to these measurements in Section 3. 

One may be puzzled as to why do we see such a stretching of space-time. Indeed, if all spatial 
distances are scaled with a universal scale factor, our local measuring units (our rulers) should also be 
stretched, and therefore we should not see the difference when comparing the two distances (e.g. the two 
wavelengths) at different times. The reason we see the difference is because we live in a gravitationally 
bound system, decoupled from the expansion of the universe: local spatial units in these systems are not 
stretched by the expansion.^ The wavelengths of photons are stretched along their geodesic path from 
one galaxy to another. In this consistent world picture, galaxies are like point particles, moving as a fluid 
in an expanding universe. 

2.1.2 The matter and energy content of the universe 

So far I have only discussed the geometrical aspects of space-time. Let us now consider the matter 
and energy content of such a universe. The most general matter fluid consistent with the assumption of 

^The subscript L refers to Luminosity, which characterizes the amount of light emitted by an object. See Eq. (|6l|). 
''The local space-time of a gravitationally bound system is described by the Schwarzschild metric, which is static [^]. 



homogeneity and isotropy is a perfect fluid, one in which an observer comoving with the fluid would see 
the universe around it as isotropic. The energy momentum tensor associated with such a fluid can be 
written as [^] 

T'^-=P5M- + (p + p);7^C/-, (7) 

where p{t) and p{t) are the pressure and energy density of the fluid at a given time in the expansion, and 
is the comoving four- velocity, satisfying U^U^ = —1. 

Let us now write the equations of motion of such a fluid in an expanding universe. According to 
general relativity, these equations can be deduced from the Einstein equations (|l|), where we substitute 
the FRW metric and the perfect fluid tensor (^. The = u = component of the Einstein equations 
constitutes the so-called Friedmann equation 
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where I have treated the cosmological constant A as a different component from matter. In fact, it can 
be associated with the vacuum energy of quantum field theory, although we still do not understand why 
should it have such a small value (120 orders of magnitude below that predicted by quantum theory), 
if it is non-zero. This constitutes today one of the most fundamental problems of physics, let alone 
cosmology. 

The conservation of energy (T^^^^, = 0), a direct consequence of the general covariance of the 
theory (G^^.^ = 0), can be written in terms of the FRW metric and the perfect fluid tensor (0) as 
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+ P^(a^)=0, (9) 



where the energy density and pressure can be split into its matter and radiation components, p = Pm+Pr, 
P = Pm + PR, with corresponding equations of state, pM = 0, pr = Pr/3. Together, the Friedmann 
and the energy-conservation equation give the evolution equation for the scale factor, 

a 47rG , , A 
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I will now make a few useful definitions. We can write the Hubble parameter today Hq in units of 
100 kms~^Mpc~^, in terms of which one can estimate the order of magnitude for the present size and 
age of the universe, 

Ho = lOO/i kms^^Mpc^S (11) 
ciJ-i = 3000/1^1 Mpc, (12) 
iJ-i = 9.773 /i^^ Gyr . (13) 

The parameter h has been measured to be in the range 0.4 < /i < 1 for decades, and only in the last few 
years has it been found to lie within 10% of h = 0.65. I will discuss those recent measurements in the 
next Section. 

One can also define a critical density pc, that which in the absence of a cosmological constant 
would correspond to a flat universe, 

Pc^^ = 1.88 10-29 g/cm^ (14) 
= 2.77 h-^ 10" MQ/{h-^ Mpc)3 , (15) 



where Mq = 1.989 x 10^^ g is a solar mass unit. The critical density p^ corresponds to approximately 
4 protons per cubic meter, certainly a very dilute fluid! In terms of the critical density it is possible to 



define the ratios Qi = Pi/pc, for matter, radiation, cosmological constant and even curvature, today. 
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We can evaluate today the radiation component JIr, corresponding to relativistic particles, from 
the density of microwave background photons, /Ocmb — fs (^'^cmb)^/(^'^)^ — 4.5 x 10~^^ g/cm'^, which 
gives f2(3MB ~ ^ 10^^ Three massless neutrinos contribute an even smaller amount. Therefore, 
we can safely neglect the contribution of relativistic particles to the total density of the universe today, 
which is dominated either by non-relativistic particles (baryons, dark matter or massive neutrinos) or by 
a cosmological constant, and write the rate of expansion in terms of its value today, 

H\a) = Hi f 4 + 4 + + 4) ■ (18) 

An interesting consequence of these redefinitions is that I can now write the Friedmann equation today, 
o = ao, as a cosmic sum rule, 

\ = Vt^^VtK^VtK, (19) 

where we have neglected f^R today. That is, in the context of a FRW universe, the total fraction of 
matter density, cosmological constant and spatial curvature today must add up to one. For instance, if 
we measure one of the three components, say the spatial curvature, we can deduce the sum of the other 
two. Making use of the cosmic sum rule today, we can write the matter and cosmological constant as a 
function of the scale factor (oq = 1) 

SvrGpM ^ / 5^ 1 

^^""^ " 3i/2(a) - a + 17M(l-a) + l^A(a3-a) \ -^-Z!- q ' ^ ^ 
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This implies that for sufficiently early times, a <^ 1, all matter-dominated FRW universes can be de- 
scribed by Einstein-de Sitter (EdS) models (i7x = 0, I^a = 0)]^ On the other hand, the vacuum energy 
will always dominate in the future. 

Another relationship which becomes very useful is that of the cosmological deceleration parameter 
today, qo, in terms of the matter and cosmological constant components of the universe, see Eq. (|Io|), 



= -^u — f^A , (22) 

^ 



which is independent of the spatial curvature. Uniform expansion corresponds to gg = and requires a 
precise cancellation: JIm = SJIa- It represents spatial sections that are expanding at a fixed rate, its scale 
factor growing by the same amount in equally-spaced time intervals. Accelerated expansion corresponds 
to go < and comes about whenever $7m < "i^k'- spatial sections expand at an increasing rate, their 
scale factor growing at a greater speed with each time interval. Decelerated expansion corresponds to 
go > and occurs whenever JIm > spatial sections expand at a decreasing rate, their scale factor 
growing at a smaller speed with each time interval. 

'Note that in the limit a ^ the radiation component starts dominating, see Eq. dla), but we still recover the EdS model. 



2.1.3 Mechanical analogy 

It is enlightening to work with a mechanical analogy of the Friedmann equation. Let us rewrite Eq. ^ 
as 

1 2 GM A 2 K 

-a a = = constant , (23) 
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where M = ^ p a is the equivalent of mass for the whole volume of the universe. Equation ( [23[ ) can 
be understood as the energy conservation law E = T + V for Si test particle of unit mass in the central 
potential 

V{r)=-^ + \kr\ (24) 
r 2 

corresponding to a Newtonian potential plus a harmonic oscillator potential with a negative spring con- 
stant k = —A/3. Note that, in the absence of a cosmological constant (A = 0), a critical universe, 
defined as the borderline between indefinite expansion and recollapse, corresponds, through the Fried- 
mann equations of motion, precisely with a flat universe {K = 0). In that case, and only in that case, a 
spatially open universe {K = — 1) corresponds to an eternally expanding universe, and a spatially closed 
universe {K = +1) to a recoUapsing universe in the future. Such a well known (textbook) correspon- 
dence is incorrect when 7^ 0: spatially open universes may recollapse while closed universes can 
expand forever. One can see in Fig. |i| a range of possible evolutions of the scale factor, for various pairs 
of values of (JIm, ^^a)- 

One can show that, for 7^ 0, a critical universe {H = H = 0) corresponds to those points 
X = oo/a > 0, for which f{x) = H'^{a) and f'{x) vanish, while f"{x) > 0, 

f{x) = x^^}M + x^nK + ^A = 0, (25) 
/'(.) = 3x^!)M + 2.n^ = { ^ : "_2si^./3si^, > „ , (26) 

Using the cosmic sum rule (p^, we can write the solutions as 

f^M < 1 
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> 1 



(28) 



The first solution corresponds to the critical point x = (a = oo), and Qk > 0, while the second one 
to X = 2|0;^|/3Qm5 and < 0. Expanding around = 1> we find JIa — 4f ~ 1)^/^M' f^i" 
r^M ^ 1- These critical solutions are asymptotic to the Einstein-de Sitter model (Qyi = 1, $7a = 0), see 
Fig. 2. 



2.1.4 Thermodynamical analogy 

It is also enlightening to find an analogy between the energy conservation equation (^ and the second 
law of Thermodynamics, 

TdS = dU + pdV , (29) 

where U = pV is the total energy of the closed system and y = is its physical volume. Equation (^ 
implies that the expansion of the universe is adiabatic or isoentropic {dS = 0), corresponding to a fluid 
in thermal equilibrium at a temperature T. For a barotropic fluid, satisfying the equation of state p = cop, 
we can write the energy density evolution as 



-{pa') = -p-{a'') = -3Hu; {pa') 



(30) 
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Fig. 1: Evolution of the scale parameter with respect to time for different values of matter density and cosmological parameter. 
The horizontal axis represents r = HQ(t — to), while the vertical axis is y — a/ao in each case. The values of (f^M, ^^a) 
for different plots are: A=(1,0), B=(0.1,0), C=(1.5,0), D=(3,0), E=(0.1,0.9), F=(0,1), G=(3,.l), H=(3,l), I=(.l,.5), J=(.5,-l), 
K=( 1-1,2.707), L=(l,2.59), M=(0.1,1.5), N=(0.1,2.5). FromRef. [ItI]. 



For relativistic particles in thermal equilibrium, the trace of the energy-momentum tensor vanishes (be- 
cause of conformal invariance) and thus Pr = Pr/3 =^ u = 1/3. In that case, the energy density of 
radiation in thermal equihbrium can be written as ^ 



PR 



i=bosons 



j=fcrmions 
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where gi^, is the number of relativistic degrees of freedom, coming from both bosons and fermions. Using 
the equilibrium expressions for the pressure and density, we can write dp = {p + p)dT/T, and therefore 



dS=^d[{p + p)V]-{p + p)v'^ 



d 



{p+p)y 

T 



const. 



(33) 



That is, up to an additive constant, the entropy per comoving volume is 5" = a^(/9 + p)V/T, which is 
conserved. The entropy per comoving volume is dominated by the contribution of relativistic particles. 
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Fig. 2: Parameter space (JIm, ^^a)- The line f^A = 1 — fiivi corresponds to a flat universe, ~ 0, separating open from 
closed universes. The line Qa = S7m/2 corresponds to uniform expansion, go ~ 0, separating accelerating from decelerating 
universes. The dashed line corresponds to critical universes, separating eternal expansion from recollapse in the future. Finally, 
the dotted line corresponds to toHo — oo, beyond which the universe has a bounce. 



SO that, to very good approximation, 

27r2 



S = -^g^siO'T)^ = constant, (34) 
g*s = J2 9t (j;) + ^ 9i (j;) ' ^^^^ 

i=bosons «=fermions 



A consequence of Eq. (34) is that, during the adiabatic expansion of the universe, the scale factor grows 
inversely proportional to the temperature of the universe, a oc 1/T. Therefore, the observational fact 
that the universe is expanding today implies that in the past the universe must have been much hotter and 
denser, and that in the future it will become much colder and dilute. Since the ratio of scale factors can 
be described in terms of the redshift parameter z, see Eq. (^, we can find the temperature of the universe 
at an earlier epoch by 

T = To (1 + z) . (36) 

Such a relation has been spectacularly confirmed with observations of absorption spectra from quasars 
at large distances, which showed that, indeed, the temperature of the radiation background scaled with 
redshift in the way predicted by the hot Big Bang model. 



2.2 Brief thermal history of the universe 

In this Section, I will briefly summarize the thermal history of the universe, from the Planck era to the 
present. As we go back in time, the universe becomes hotter and hotter and thus the amount of energy 
available for particle interactions increases. As a consequence, the nature of interactions goes from those 
described at low energy by long range gravitational and electromagnetic physics, to atomic physics, 
nuclear physics, all the way to high energy physics at the electroweak scale, gran unification (perhaps), 
and finally quantum gravity. The last two are still uncertain since we do not have any experimental 
evidence for those ultra high energy phenomena, and perhaps Nature has followed a different path. f\ 

*See the recent theoretical developments on large extra dimensions and quantum gravity at the TeV \w. 



The way we know about the high energy interactions of matter is via particle accelerators, which 
are unravelling the details of those fundamental interactions as we increase in energy. However, one 
should bear in mind that the physical conditions that take place in our high energy colliders are very 
different from those that occurred in the early universe. These machines could never reproduce the 
conditions of density and pressure in the rapidly expanding thermal plasma of the early universe. Nev- 
ertheless, those experiments are crucial in understanding the nature and rate of the local fundamental 
interactions available at those energies. What interests cosmologists is the statistical and thermal proper- 
ties that such a plasma should have, and the role that causal horizons play in the final outcome of the early 
universe expansion. For instance, of crucial importance is the time at which certain particles decoupled 
from the plasma, i.e. when their interactions were not quick enough compared with the expansion of the 
universe, and they were left out of equilibrium with the plasma. 

One can trace the evolution of the universe from its origin till today. There is still some speculation 
about the physics that took place in the universe above the energy scales probed by present colliders. 
Nevertheless, the overall layout presented here is a plausible and hopefully testable proposal. According 
to the best accepted view, the universe must have originated at the Planck era (10^^ GeV, 10~^^ s) 
from a quantum gravity fluctuation. Needless to say, we don't have any experimental evidence for such 
a statement: Quantum gravity phenomena are still in the realm of physical speculation. However, it 
is plausible that a primordial era of cosmological inflation originated then. Its consequences will be 
discussed below. Soon after, the universe may have reached the Grand Unified Theories (GUT) era (10^^ 
GeV, 10~^^ s). Quantum fluctuations of the inflaton field most probably left their imprint then as tiny 
perturbations in an otherwise very homogenous patch of the universe. At the end of inflation, the huge 
energy density of the inflaton field was converted into particles, which soon thermalized and became the 
origin of the hot Big Bang as we know it. Such a process is called reheating of the universe. Since 
then, the universe became radiation dominated. It is probable (although by no means certain) that the 
asymmetry between matter and antimatter originated at the same time as the rest of the energy of the 
universe, from the decay of the inflaton. This process is known under the name of baryogenesis since 
baryons (mostly quarks at that time) must have originated then, from the leftovers of their annihilation 
with antibaryons. It is a matter of speculation whether baryogenesis could have occurred at energies 
as low as the electroweak scale (100 GeV, 10"^'' s). Note that although particle physics experiments 
have reached energies as high as 100 GeV, we still do not have observational evidence that the universe 
actually went through the EW phase transition. If confirmed, baryogenesis would constitute another 
"window" into the early universe. As the universe cooled down, it may have gone through the quark- 
gluon phase transition (10^ MeV, 10~^ s), when baryons (mainly protons and neutrons) formed from 
their constituent quarks. 

The furthest window we have on the early universe at the moment is that of primordial nucleosyn- 
thesis (1 — 0.1 MeV, I s - 3 min), when protons and neutrons were cold enough that bound systems 
could form, giving rise to the lightest elements, soon after neutrino decoupling: It is the realm of nuclear 
physics. The observed relative abundances of light elements are in agreement with the predictions of the 
hot Big Bang theory. Immediately afterwards, electron-positron annihilation occurs (0.5 MeV, 1 min) 
and all their energy goes into photons. Much later, at about (I eV, ~ 10^ yr), matter and radiation have 
equal energy densities. Soon after, electrons become bound to nuclei to form atoms (0.3 eV, 3 x 10^ 
yr), in a process known as recombination: It is the realm of atomic physics. Immediately after, photons 
decouple from the plasma, travelling freely since then. Those are the photons we observe as the cosmic 
microwave background. Much later (~ 1 — 10 Gyr), the small inhomogeneities generated during inflation 
have grown, via gravitational collapse, to become galaxies, clusters of galaxies, and superclusters, char- 
acterizing the epoch of structure formation. It is the realm of long range gravitational physics, perhaps 
dominated by a vacuum energy in the form of a cosmological constant. Finally (3K, 13 Gyr), the Sun, 
the Earth, and biological Ufe originated from previous generations of stars, and from a primordial soup 
of organic compounds, respectively. 



I will now review some of the more robust features of the Hot Big Bang theory of which we have 
precise observational evidence. 

2.2.1 Primordial nucleosynthesis and light element abundance 

In this subsection I will briefly review Big Bang nucleosynthesis and give the present observational 
constraints on the amount of baryons in the universe. In 1920 Eddington suggested that the sun might 
derive its energy from the fusion of hydrogen into helium. The detailed reactions by which stars burn 
hydrogen were first laid out by Hans Bethe in 1939. Soon afterwards, in 1946, George Gamow realized 
that similar processes might have occurred also in the hot and dense early universe and gave rise to the 
first light elements iQ]. These processes could take place when the universe had a temperature of around 
Tj^g ~ 1 — 0.1 MeV, which is about 100 times the temperature in the core of the Sun, while the density 
is /?^3 = fo5*^Ns ~ gcm~^, about the same density as the core of the Sun. Note, however, that 
although both processes are driven by identical thermonuclear reactions, the physical conditions in star 
and Big Bang nucleosynthesis are very different. In the former, gravitational collapse heats up the core of 
the star and reactions last for billions of years (except in supernova explosions, which last a few minutes 
and creates all the heavier elements beyond iron), while in the latter the universe expansion cools the hot 
and dense plasma in just a few minutes. Nevertheless, Gamow reasoned that, although the early period of 
cosmic expansion was much shorter than the lifetime of a star, there was a large number of free neutrons 
at that time, so that the lighter elements could be built up quickly by succesive neutron captures, starting 
with the reaction n + p ^ D + ^. The abundances of the light elements would then be correlated with 
their neutron capture cross sections, in rough agreement with observations [T0[ ]. 
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Fig. 3: The relative abundance of light elements to Hidrogen. Note the large range of scales involved. From Ref. m: 



Nowadays, Big Bang nucleosynthesis (BBN) codes compute a chain of around 30 coupled nuclear 



reactions, to produce all the light elements up to beryllium-7. f\ Only the first four or five elements can 
be computed with accuracy better than 1% and compared with cosmological observations. These light 
elements are H, ^He, D, ^He, '^Li, and perhaps also ^Li. Their observed relative abundance to hydrogen 
is [1 : 0.25 : 3 • 10"^ : 2 • 10~^ : 2 • 10"^°] with various errors, mainly systematic. The BBN 
codes calculate these abundances using the laboratory measured nuclear reaction rates, the decay rate of 
the neutron, the number of light neutrinos and the homogeneous FRW expansion of the universe, as a 
function of only one variable, the number density fraction of baryons to photons, r] = UB/rij. In fact, 
the present observations are only consistent, see Fig. I and Ref. ^ with a very narrow range of 
values of 

7710 = 10^% = 4.6 - 5.9 . (37) 

Such a small value of t] indicates that there is about one baryon per 10^ photons in the universe today. 
Any acceptable theory of baryogenesis should account for such a small number. Furthermore, the present 
baryon fraction of the critical density can be calculated from rjiQ as [ |To| ] 

nBh^ = 3.6271 X 10"^ r/io = 0.0190 ± 0.0024 (95% c.l.) (38) 

Clearly, this number is well below closure density, so baryons cannot account for all the matter in the 
universe, as I shall discuss below. 



2.2.2 Neutrino decoupling 

Just before the nucleosynthesis of the lightest elements in the early universe, weak interactions were too 
slow to keep neutrinos in thermal equilibrium with the plasma, so they decoupled. We can estimate the 
temperature at which decoupling occurred from the weak interaction cross section, a„ ~ CpT"^ at finite 
temperature T, where Gp = 1.2 x 10~^ GeV^^ is the Fermi constant. The neutrino interaction rate, via 
W boson exchange in n + u ^ p + and p + i> ^ n + e~^, can be written as [||] 

r, = n,{a^\v\) ^Gj,T\ (39) 

while the rate of expansion of the universe at that time (g^: = 10.75) was H ~ 5.4 T^/Mp, where 
Mp = 1.22 X 10^^ GeV is the Planck mass. Neutrinos decouple when their interaction rate is slower 
than the universe expansion, Ty < H or, equivalently, at Ty^^cc — 0.8 MeV. Below this temperature, 
neutrinos are no longer in thermal equilibrium with the rest of the plasma, and their temperature continues 
to decay inversely proportional to the scale factor of the universe. Since neutrinos decoupled before 
e+e^ annihilation, the cosmic background of neutrinos has a temperature today lower than that of the 
microwave background of photons. Let us compute the difference. At temperatures above the the mass 
of the electron, T > nif. = 0.511 MeV, and below 0.8 MeV, the only particle species contributing to 
the entropy of the universe are the photons {g^ = 2) and the electron-positron pairs {g^ = 4 x |); total 
number of degrees of freedom (7* = ^. At temperatures T ~ mg, electrons and positrons annihilate into 
photons, heating up the plasma (but not the neutrinos, which had decoupled already). At temperatures 
T < rrie, only photons contribute to the entropy of the universe, with g^, = 2 degrees of freedom. 
Therefore, from the conservation of entropy, we find that the ratio of Ty and Ty today must be 

^ = (t)'^'"^-^°^ ^ = 1.945 K, (40) 

where I have used T^^g = 2.725 ± 0.002 K. We still have not measured such a relic background of 
neutrinos, and probably will remain undetected for a long time, since they have an average energy of 
order 10^^ eV, much below that required for detection by present experiments (of order GeV), precisely 
because of the relative weakness of the weak interactions. Nevertheless, it would be fascinating if, in the 
future, ingenious experiments were devised to detect such a background, since it would confirm one of 
the most robust features of Big Bang cosmology. 

'The rest of nuclei, up to iron (Fe), are produced in lieavy stars, and beyond Fe in novae and supemovae explosions. 



2.2.3 Matter-radiation equality 

Relativistic species have energy densities proportional to the quartic power of temperature and therefore 
scale as /Or oc a~^, while non-relativistic particles have essentially zero pressure and scale as py[ oc a"^, 



see Eq. (30). Therefore, there will be a time in the evolution of the universe in which both energy 
densities are equal /7R(teq) = /OM(ieq)- Since then both decay differently, and thus 
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3.1 X 10^ Om/i^ 



(41) 



where I have used I^r/i^ = Q.^^^h? + Q.yh^ = 3.24 x 10 ^ for three massless neutrinos at T = Ty. As 
I will show later, the matter content of the universe today is below critical, JIm — 0.3, while h ~ 0.65, 
and therefore (1 + Zgq) ~ 3900, or about icq = 1-2 x 10^ {Q.Mh'^y^ ~ 7 x 10^ years after the origin of 



the universe. Around the time of matter-radiation equality, the rate of expansion ( |18[ ) can be written as 

(ao = 1) 



H{a) = Ho ( a ^ + I^m a 
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(42) 



The horizon size is the coordinate distance travelled by a photon since the beginning of the universe, 
dn ~ H^^, i.e. the size of causally connected regions in the universe. The comoving horizon size is 
then given by 

aH[a) u ivi y ^ y 

Thus the horizon size at matter-radiation equality (a = Ocq) is 



(44) 



This scale plays a very important role in theories of structure formation. 




Fig. 4: The equilibrium ionization fraction X^^' as a function of redstiift. Tiie two lines show the range of 7710 — 4.6 — 5.9. 



2.2.4 Recombination and photon decoupling 

As the temperature of the universe decreased, electrons could eventually become bound to protons to 
form neutral hydrogen. Nevertheless, there is always a non-zero probability that a rare energetic photon 
ionizes hydrogen and produces a free electron. The ionization fraction of electrons in equilibrium with 
the plasma at a given temperature is given by |^ 

1-^1" _ 4^/2C(3) ^ ^^^3/2 ^^^^^^^ 



rj I — e^-°/' , (45) 



where ii^ion = 13.6 eV is the ionization energy of hydrogen, and rj is the baryon-to-photon ratio (|37|). If 
we now use Eq. (^), we can compute the ionization fraction X^"^ as a function of redshift z, see Fig. ^ 
Note that the huge number of photons with respect to electrons (in the ratio ^He : : 7 ~ 1 : 4 : 10^'^) 
implies that even at a very low temperature, the photon distribution will contain a sufficiently large num- 
ber of high-energy photons to ionize a significant fraction of hydrogen. In fact, defining recombination 
as the time at which X^'^ = 0.1, one finds that the recombination temperature is Tree = 0.3 eV <C ^^ion> 
for r/10 ^ 5.2. Comparing with the present temperature of the microwave background, we deduce the 
corresponding redshift at recombination, (1 + Zrec) ~ 1270. 

Photons remain in thermal equilibrium with the plasma of baryons and electrons through elastic 
Thomson scattering, with cross section 



3m^ 



6.65 X 10"^^ cm^ = 0.665 barn, 



(46) 



where a = 1/137.036 is the dimensionless electromagnetic coupling constant. The mean free path of 
photons in such a plasma can be estimated from the photon interaction rate, ~ = nefjj,. 
For temperatures above a few eV, the mean free path is much smaller that the causal horizon at that 
time and photons suffer multiple scattering: the plasma is like a dense fog. Photons will decouple from 
the plasma when their interaction rate cannot keep up with the expansion of the universe and the mean 
free path becomes larger than the horizon size: the universe becomes transparent. We can estimate 
this moment by evaluating = H at photon decoupling. Using Ue = X(.r] riy, one can compute 
the decoupling temperature as Tdoc = 0.26 eV, and the corresponding redshift as (1 + Zdec) — 1100. 
This redshift defines the so called last scattering surface, when photons last scattered off protons and 
electrons and travelled freely ever since. This decoupling occurred when the universe was approximately 
tdec = 1-8 X 10^ (17m^^)"^/^ ~ 5 X 10^ years old. 
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Fig. 5: The Cosmic Microwave Background Spectrum seen by the HRAS instrument on COBE. The left panel corresponds to 
the monopole spectrum, To = 2.725 ± 0.002 K, where the error bars are smaller than the line width. The right panel shows the 
dipole spectrum, 6T1 = 3.372 ± 0.014 mK. From Ref. [Ol. 



2.2.5 The microwave background 

One of the most remarkable observations ever made my mankind is the detection of the relic background 
of photons from the Big Bang. This background was predicted by George Gamow and collaborators 
in the 1940s, based on the consistency of primordial nucleosynthesis with the observed helium abun- 
dance. They estimated a value of about 10 K, although a somewhat more detailed analysis by Alpher and 
Herman in 1950 predicted Tj ^ 5K. Unfortunately, they had doubts whether the radiation would have 



survived until the present, and this remarkable prediction slipped into obscurity, until Dicke, Peebles, 



Roll and Wilkinson [ ]13| ] studied the problem again in 1965. Before they could measure the photon back- 
ground, they learned that Penzias and Wilson had observed a weak isotropic background signal at a radio 
wavelength of 7.35 cm, corresponding to a blackbody temperature of = 3.5 it 1 K. They published 
their two papers back to back, with that of Dicke et al. explaining the fundamental significance of their 
measurement 

Since then many different experiments have confirmed the existence of the microwave background. 
The most outstanding one has been the Cosmic Background Explorer (COBE) satellite, whose FIRAS 
instrument measured the photon background with great accuracy over a wide range of frequencies {u = 
1 — 97 cm^^), see Ref. [[TJ], with a spectral resolution — = 0.0035. Nowadays, the photon spectrum is 
confirmed to be a blackbody spectrum with a temperature given by [p^ 



- CMB 



2.725 ± 0.002 K (systematic, 95% c.l.) ± 7 /iK (Icr statistical) 



(47) 



In fact, this is the best blackbody spectrum ever measured, see Fig. ^, with spectral distortions below the 
level of 10 parts per million (ppm). 
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Fig. 6: The Cosmic Microwave Background Spectrum seen by the DMR instrument on COBE. The top figure corresponds to 
the monopole, To = 2.725 ± 0.002 K. The middle figure shows the dipole, 6Ti = 3.372 ± 0.014 mK, and the lower figure 
shows the quadrupole and higher multipoles, 5T2 = 18 ± 2 ^K. The central region corresponds to foreground by the galaxy. 
From Ref. 10. 

Moreover, the differential microwave radiometer (DMR) instrument on COBE, with a resolution 
of about 7° in the sky, has also confirmed that it is an extraordinarily isotropic background. The devia- 
tions from isotropy, i.e. differences in the temperature of the blackbody spectrum measured in different 



directions in the sky, are of the order of 20 /xK on large scales, or one part in 10^, see Ref. [14]. There 
is, in fact, a dipole anisotropy of one part in 10^, 5Ti = 3.372 it 0.007 mK (95% c.l.), in the direction 
of the Virgo cluster, (/, b) = (264.14° ± 0.30, 48.26° ± 0.30) (95% c.l.). Under the assumption that a 
Doppler effect is responsible for the entire CMB dipole, the velocity of the Sun with respect to the CMB 



rest frame is vq = 371 it 0.5 km/s, see Ref. [12].0 When subtracted, we are left with a whole spectrum 
of anisotropics in the higher multipoles (quadrupole, octupole, etc.), ST2 = 18 it 2 /xK (95% c.L), see 



Ref. [|14p and Fig. g. 

Soon after COBE, other groups quickly confirmed the detection of temperature anisotropics at 
around 30 /xK and above, at higher multipole numbers or smaller angular scales. As I shall discuss below, 
these anisotropics play a crucial role in the understanding of the origin of structure in the universe. 

2.3 Large-scale structure formation 

Although the isotropic microwave background indicates that the universe in the past was extraordinarily 
homogeneous, we know that the universe today is not exactly homogeneous: we observe galaxies, clus- 
ters and superclusters on large scales. These structures are expected to arise from very small primordial 
inhomogeneities that grow in time via gravitational instability, and that may have originated from tiny 
ripples in the metric, as matter fell into their troughs. Those ripples must have left some trace as temper- 
ature anisotropics in the microwave background, and indeed such anisotropics were finally discovered 
by the COBE satellite in 1992. The reason why they took so long to be discovered was that they appear 
as perturbations in temperature of only one part in 10^. 

While the predicted anisotropics have finally been seen in the CMB, not all kinds of matter and/or 
evolution of the universe can give rise to the structure we observe today. If we define the density contrast 



as 1 15] 



5(f , a) = ^(^'«) /^(q) = / d-'k 6k{a) e'^-"-' , (48) 



p{a) 

where p{a) = po a^'^ is the average cosmic density, we need a theory that will grow a density contrast 
with amplitude 6 ~ 10^^ at the last scattering surface (z = 1100) up to density contrasts of the order of 
5 ~ 10^ for galaxies at redshifts z <^ 1, i.e. today. This is a necessary requirement for any consistent 



theory of structure formation [16|. 

Furthermore, the anisotropics observed by the COBE satellite correspond to a small-amplitude 
scale-invariant primordial power spectrum of inhomogeneities 

P{k) = (|(5fcp) oc A:" , with n = 1 , (49) 

where the brackets (•) represent integration over an ensemble of different universe realizations. These 
inhomogeneities are like waves in the space-time metric. When matter fell in the troughs of those waves, 
it created density perturbations that collapsed gravitationally to form galaxies and clusters of galaxies, 
with a spectrum that is also scale invariant. Such a type of spectrum was proposed in the early 1970s by 
Edward R. Harrison, and independently by the Russian cosmologist Yakov B. Zel'dovich, see Ref. []T7[], 
to explain the distribution of galaxies and clusters of galaxies on very large scales in our observable 
universe. 

Today various telescopes - Uke the Hubble Space Telescope, the twin Keck telescopes in Hawaii 
and the European Southern Observatory telescopes in Chile - are exploring the most distant regions 
of the universe and discovering the first galaxies at large distances. The furthest galaxies observed so 
far are at redshifts of z ~ 5, or 12 billion light years from the Earth, whose light was emitted when 
the universe had only about 5% of its present age. Only a few galaxies are known at those redshifts, 
but there are at present various catalogs like the CfA and APM galaxy catalogs, and more recently the 
IRAS Point Source redshift Catalog, see Fig. 0, and Las Campanas redshift surveys, that study the spatial 
distribution of hundreds of thousands of galaxies up to distances of a billion light years, or z < 0.1, that 
recede from us at speeds of tens of thousands of kilometres per second. These catalogs are telling us 
about the evolution of clusters of galaxies in the universe, and already put constraints on the theory of 



COBE even determined the annual variation due to the Earth's motion around the Sun - the ultimate proof of Copernicus' 
hypothesis. 
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Fig. 7: The IRAS Point Source Catalog redshift survey contains some 15,000 galaxies, covering over 83% of the sky up to 
redshifts of z < 0.05. We show here the projection of the galaxy distribution in galactic coordinates. From Ref. [|l8||. 



structure formation. From these observations one can infer that most galaxies formed at redshifts of the 
order of 2 — 6; clusters of galaxies formed at redshifts of order 1, and superclusters are forming now. 
That is, cosmic structure formed from the bottom up: from galaxies to clusters to superclusters, and not 
the other way around. This fundamental difference is an indication of the type of matter that gave rise 
to structure. The observed power spectrum of the galaxy matter distribution from a selection of deep 
redshift catalogs can be seen in Fig. ^. 

We know from Big Bang nucleosynthesis that all the baryons in the universe cannot account for 
the observed amount of matter, so there must be some extra matter (dark since we don't see it) to account 
for its gravitational pull. Whether it is relativistic (hot) or non-relativistic (cold) could be inferred from 
observations: relativistic particles tend to diffuse from one concentration of matter to another, thus trans- 
ferring energy among them and preventing the growth of structure on small scales. This is excluded by 
observations, so we conclude that most of the matter responsible for structure formation must be cold. 
How much there is is a matter of debate at the moment. Some recent analyses suggest that there is not 
enough cold dark matter to reach the critical density required to make the universe flat. If we want to 
make sense of the present observations, we must conclude that some other form of energy permeates 
the universe. In order to resolve this issue, even deeper galaxy redshift catalogs are underway, looking 
at millions of galaxies, like the Sloan Digital Sky Survey (SDSS) and the Anglo-Australian two degree 
field (2dF) Galaxy Redshift Survey, which are at this moment taking data, up to redshifts of z ^ 0.5, 
over a large region of the sky. These important observations will help astronomers determine the nature 
of the dark matter and test the validity of the models of structure formation. 

Before COBE discovered the anisotropics of the microwave background there were serious doubts 
whether gravity alone could be responsible for the formation of the structure we observe in the universe 
today. It seemed that a new force was required to do the job. Fortunately, the anisotropics were found 
with the right amplitude for structure to be accounted for by gravitational collapse of primordial inho- 
mogeneities under the attraction of a large component of non-relativistic dark matter. Nowadays, the 
standard theory of structure formation is a cold dark matter model with a non vanishing cosmological 
constant in a spatially flat universe. Gravitational collapse amplifies the density contrast initially through 
linear growth and later on via non-linear collapse. In the process, overdense regions decouple from 
the Hubble expansion to become bound systems, which start attracting eachother to form larger bound 
structures. In fact, the largest structures, superclusters, have not yet gone non-linear. 

The primordial spectrum (49) is reprocessed by gravitational instability after the universe becomes 
matter dominated and inhomogeneities can grow. Linear perturbation theory shows that the growing 
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Fig. 8: The left panel shows the matter power spectrum for clusters of galaxies, from three different cluster surveys. The right 
panel shows a compilation of the most recent estimates of the power spectrum of galaxy clustering, from four of the largest 
available redshift surveys of optically-selected galaxies, compared to the deprojected spectrum of the 2D APM galaxy survey. 
FromRef. 



mode of small density contrasts go like [15,16] 



6{a) oc a 
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a > a, 



eq 
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(50) 



in the Einstein-de Sitter limit (lo = p/ p = 1/3 and 0, for radiation and matter, respectively). There are 
slight deviations for a » Ceq, if / 1 or J7a / 0, but we will not be concerned with them here. 
The important observation is that, since the density contrast at last scattering is of order 6 10^^, and 
the scale factor has grown since then only a factor zjec ~ 10^, one would expect a density contrast 
today of order ~ 10"^. Instead, we observe structures like galaxies, where 6 ~ 10^. So how can 
this be possible? The microwave background shows anisotropics due to fluctuations in the baryonic 
matter component only (to which photons couple, electromagnetically). If there is an additional matter 
component that only couples through very weak interactions, fluctuations in that component could grow 
as soon as it decoupled from the plasma, well before photons decoupled from baryons. The reason why 
baryonic inhomogeneities cannot grow is because of photon pressure: as baryons collapse towards denser 
regions, radiation pressure eventually halts the contraction and sets up acoustic oscillations in the plasma 
that prevent the growth of perturbations, until photon decoupling. On the other hand, a weakly interacting 
cold dark matter component could start gravitational collapse much earlier, even before matter-radiation 
equality, and thus reach the density contrast amplitudes observed today. The resolution of this mismatch 
is one of the strongest arguments for the existence of a weakly interacting cold dark matter component 
of the universe. 

How much dark matter there is in the universe can be deduced from the actual power spectrum (the 
Fourier transform of the two-point correlation function of density perturbations) of the observed large 
scale structure. One can decompose the density contrast in Fourier components, see Eq. (|4^. This is 
very convenient since in linear perturbation theory individual Fourier components evolve independently. 
A comoving wavenumber k is said to "enter the horizon" when k = djj^{a) = aH{a). If a certain 
perturbation, of wavelength A = k^^ < dH{acq), enters the horizon before matter-radiation equality, the 
fast radiation-driven expansion prevents dark-matter perturbations from collapsing. Since light can only 
cross regions that are smaller than the horizon, the suppression of growth due to radiation is restricted 
to scales smaller than the horizon, while large-scale perturbations remain unaffected. This is the reason 
why the horizon size at equality, Eq. (p4|), sets an important scale for structure growth. 



/cgq = d^^(aeq) ^ 0.083 (nuh) h Mpc 
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(51) 



'The decaying modes go like 5{t) t ^, for all uj. 
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Fig. 9: The power spectrum for cold dark matter (CDM), tilted cold dark matter (TCDM), hot dark matter (HDM), and mixed 
hot plus cold dark matter (MDM), normalized to COBE, for large-scale structure formation. From Ref. |^o|. 



The suppression factor can be easily computed from ( |50| ) as /sup 
words, the processed power spectrum P{k) will have the form: 



(aenter/Oeq)^ = {kcq/kf- In Other 
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This is precisely the shape that large-scale galaxy catalogs are bound to test in the near future, see Fig. ^ 
Furthermore, since relativistic Hot Dark Matter (HDM) transfer energy between clumps of matter, they 
will wipe out small scale perturbations, and this should be seen as a distinctive signature in the matter 
power spectra of future galaxy catalogs. On the other hand, non-relativistic Cold Dark Matter (CDM) 
allow structure to form on all scales via gravitational collapse. The dark matter will then pull in the 
baryons, which will later shine and thus allow us to see the galaxies. 

Naturally, when baryons start to collapse onto dark matter potential wells, they will convert a large 
fraction of their potential energy into kinetic energy of protons and electrons, ionizing the medium. As a 
consequence, we expect to see a large fraction of those baryons constituting a hot ionized gas surrounding 
large clusters of galaxies. This is indeed what is observed, and confirms the general picture of structure 
formation. 



3 DETERMINATION OF COSMOLOGICAL PARAMETERS 

In this Section, I will restrict myself to those recent measurements of the cosmological parameters by 
means of standard cosmological techniques, together with a few instances of new results from recently 
applied techniques. We will see that a large host of observations are determining the cosmological 
parameters with some reliability of the order of 10%. However, the majority of these measurements are 
dominated by large systematic errors. Most of the recent work in observational cosmology has been 
the search for virtually systematic-free observables, like those obtained from the microwave background 
anisotropics, and discussed in Section 4.4. I will devote, however, this Section to the more 'classical' 
measurements of the following cosmological parameters: The rate of expansion Hq; the matter content 
r^M; the cosmological constant JIa; the spatial curvature r2x> and the age of the universe toj3 

These five basic cosmological parameters are not mutually independent. Using the homogeneity 
and isotropy on large scales observed by COBE, we can infer relationships between the different cosmo- 

'"We will take the baryon fraction as given by observations of light element abundances, in accordance with Big Bang 
nucleosynthesis, see Eq. (138). 



logical parameters through the Einstein-Friedmann equations. In particular, we can deduce the value of 
the spatial curvature from the Cosmic Sum Rule, 



1 = + ^^A + : 



(53) 



or viceversa, if we determine that the universe is spatially flat from observations of the microwave back- 
ground, we can be sure that the sum of the matter content plus the cosmological constant must be one. 

Another relationship between parameters appears for the age of the universe. In a FRW cosmology, 
the cosmic expansion is determined by the Friedmann equation (js]). Defining a new time and normaUzed 
scale factor, 

y = — = , r = Ho{t - to) , (54) 

ao 1 + z 

we can write the Friedmann equation with the help of the Cosmic Sum Rule ( |T9| ) as 



y'{T) = \l + (y-i - 1)17m + (y' - 1)5^/ 
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(55) 



with initial condition y(0) = 1, y'(0) = 1. Therefore, the present age to is a function of the other 
parameters, to = /{Hq, ^m, ^a), determined from 
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We show in Fig. 11^ the contour hnes for constant to-f^o in parameter space (Om, ^^a)- 



(56) 
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Fig. 10: The contour lines correspond to equal toHo = 0.5 — 1.0, 1.2, 1.5, 2.0 and 5.0, from bottom to top, in parameter space 
(r^M, ^^a). The hne toHo = oo would be indistinguishable from that of toHo = 5. From Ref. |t7||. 



There are two specific limits of interest: an open universe with Q\ = 0, for which the age is given by 
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and aflat universe with JIa = 1 — ^M, for which the age can also be expressed in compact form. 
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We have plotted these functions in Fig. 11. It is clear that in both cases to^^o 2/3 as 1- 
We can now use these relations as a consistency check between the cosmological observations of Hq, 
Qm, and to. Of course, we cannot measure the age of the universe directly, but only the age of its 
constituents: stars, galaxies, globular clusters, etc. Thus we can only find a lower bound on the age of 
the universe, to ~ ^gai + 1-5 Gyr. As we will see, this is not a trivial bound and, in several occasions, 
during the progress towards better determinations of the cosmological parameters, the universe seemed 
to be younger than its constituents, a logical inconsistency, of course, only due to an incorrect assessment 



of systematic errors [21]. 
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Fig. 1 1 : The age of the universe as a function of the matter content, for an open and a flat universe. From Ref. ||22|]. 

In order to understand those recent measurements, one should also define what is known as the 
luminosity distance to an object in the universe. Imagine a source that is emitting light at a distance 
from a detector of area dA. The absolute luminosity C of such a source is nothing but the energy emitted 
per unit time. A standard candle is a luminous object that can be calibrated with some accuracy and 
therefore whose absolute luminosity is known, within certain errors. For example, Cepheid variable stars 
and type la supernovae are considered to be reasonable standard candles, i.e. their calibration errors are 
within bounds. The energy flux received at the detector is the measured energy per unit time per unit 
area of the detector coming from that source. The luminosity distance cZl is then defined as the radius 
of the sphere centered on the source for which the absolute luminosity would give the observed flux, 
= C/AiTd\. In a Friedmann-Robertson-Walker universe, light travels along null geodesies, ds^ = 0, 
or, see Eq. (Q), 

dr 1 dz ^^^^ 



1 + alHl r2 '^0^0 V(l + zf{l + z^u) - z{2 + z)Oa ' 

which determines the coordinate distance r = r{z, Hq, Q,m, ^a), as a function of redshift z and the other 
cosmological parameters. Now let us consider the effect of the universe expansion on the observed flux 
coming from a source at a certain redshift z from us. First, the photon energy on its way here will be 
redshifted, and thus the observed energy Eq = E/{\ + z). Second, the rate of photon arrival will be 
time-delayed with respect to that emitted by the source, dto = (1 + z)dt. Finally, the fraction of the area 
of the 2-sphere centered on the source that is covered by the detector is dA/A-na^ r'^{z). Therefore, the 
total flux detected is 

C _ C 

^ " 47ra2r2(z) = 4^ ' ^^^^ 



The final expression for tiie luminosity distance di, as a function of redshift is thus given by [||] 
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where sinn(j;) = x K = 0; sin(x) if i^' = +1 and sinh(x) if K 
around z = 0, we obtain Eq. 



-1. Expanding to second order 
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This expression goes beyond the leading linear term, corresponding to the Hubble law, into the second 
order term, which is sensitive to the cosmological parameters Um and Q\. It is only recently that cos- 
mological observations have gone far enough back into the early universe that we can begin to probe the 
second term, as I will discuss shortly. Higher order terms are not yet probed by cosmological observa- 
tions, but they would contribute as important consistency checks. 

Let us now pursue the analysis of the recent determinations of the most important cosmological 
parameters: the rate of expansion Hq, the matter content ilM> the cosmological constant the spatial 
curvature D.k, and the age of the universe to- 



3.1 The rate of expansion Hq 



Over most of last century the value of Hq has been a constant source of disagreement [|21[]. Around 
1929, Hubble measured the rate of expansion to be Hq = 500 kms^^Mpc^^, which implied an age of 
the universe of order to ~ 2 Gyr, in clear conflict with geology. Hubble's data was based on Cepheid 
standard candles that were incorrectly calibrated with those in the Large Magellanic Cloud. Later on, 
in 1954 Baade recalibrated the Cepheid distance and obtained a lower value, Hq = 250 kms^^Mpc^^, 
still in conflict with ratios of certain unstable isotopes. Finally, in 1958 Sandage realized that the bright- 
est stars in galaxies were ionized HII regions, and the Hubble rate dropped down to Hq = 60 kms^^ 
Mpc^^, still with large (factor of two) systematic errors. Fortunately, in the past 15 years there has 
been significant progress towards the determination of Hq, with systematic errors approaching the 10% 
level. These improvements come from two directions. First, technological, through the replacement of 
photographic plates (almost exclusively the source of data from the 1920s to 1980s) with charged couple 
devices (CCDs), i.e. solid state detectors with excellent flux sensitivity per pixel, which were previously 
used successfully in particle physics detectors. Second, by the refinement of existing methods for mea- 
suring extragalactic distances (e.g. parallax, Cepheids, supernovae, etc.). Finally, with the development 
of completely new methods to determine Hq, which fall into totally independent and very broad cate- 
gories: a) Gravitational lensing; b) Sunyaev-Zel'dovich effect; c) Extragalactic distance scale, mainly 
Cepheid variability and type la Supernovae; d) Microwave background anisotropics. I will review here 
the first three, and leave the last method for Section 4.4, since it involves knowledge about the primordial 
spectrum of inhomogeneities. 



3.1.1 Gravitational lensing 

Imagine a quasi-stellar object (QSO) at large redshift (z » 1) whose light is lensed by an intervening 
galaxy at redshift z ~ 1 and arrives to an observer at z = 0. There will be at least two different 
images of the same background variable point source. The arrival times of photons from two different 
gravitationally lensed images of the quasar depend on the different path lengths and the gravitational 
potential traversed. Therefore, a measurement of the time delay and the angular separation of the different 
images of a variable quasar can be used to determine Hq with great accuracy. This method, proposed in 



1964 by Refsdael [23], offers tremendous potential because it can be applied at great distances and it is 



based on very solid physical principles | |24| ] 



Unfortunately, there are very few systems with both a favourable geometry (i.e. a known mass 
distribution of the intervening galaxy) and a variable background source with a measurable time delay. 
That is the reason why it has taken so much time since the original proposal for the first results to come 
out. Fortunately, there are now very powerful telescopes that can be used for these purposes. The best 
candidate to-date is the QSO 0957 + 561, observed with the 10m Keck telescope, for which there is a 
model of the lensing mass distribution that is consistent with the measured velocity dispersion. Assuming 
a flat space with Qm = 0.25, one can determine pS] ] 

ifo = 72 ± 7 (Icr statistical) ± 15% (systematic) kms~^Mpc"^ . (63) 

The main source of systematic error is the degeneracy between the mass distribution of the lens and 
the value of Hq. Knowledge of the velocity dispersion within the lens as a function of position helps 
constrain the mass distribution, but those measurements are very difficult and, in the case of lensing by 
a cluster of galaxies, the dark matter distribution in those systems is usually unknown, associated with a 
complicated cluster potential. Nevertheless, the method is just starting to give promising results and, in 
the near future, with the recent discovery of several systems with optimum properties, the prospects for 
measuring Hq and lowering its uncertainty with this technique are excellent. 



3.1.2 Sunyaev-Zel'dovich effect 

As discussed in the previous Section, the gravitational collapse of baryons onto the potential wells gen- 
erated by dark matter gave rise to the reionization of the plasma, generating an X-ray halo around rich 
clusters of galaxies, see Fig. 12. The inverse-Compton scattering of microwave background photons 



off the hot electrons in the X-ray gas results in a measurable distortion of the blackbody spectrum 
of the microwave background, known as the Sunyaev-Zel'dovich (SZ) effect. Since photons acquire 
extra energy from the X-ray electrons, we expect a shift towards higher frequencies of the spectrum, 
(Av/v) ~ {kBTgas / nT-eC^) ~ 10^^. This corresponds to a Jecremen? of the microwave background tem- 
perature at low frequencies (Rayleigh- Jeans region) and an increment at high frequencies, see Ref. [26|. 





Fig. 12: The Coma cluster of galaxies, seen here in an optical image (left) and an X-ray image (right), taken by the recently 
launched Chandra X-ray Observatory. From Ref. pv]]. 

Measuring the spatial distribution of the SZ effect (3 K spectrum), together with a high resolution 
X-ray map (10* K spectrum) of the cluster, one can determine the density and temperature distribution 
of the hot gas. Since the X-ray flux is distance-dependent (J^ = £/A'Kd\), while the SZ decrement is 
not (because the energy of the CMB photons increases as we go back in redshift, v = i'o{l + z), and 
exactly compensates the redshift in energy of the photons that reach us), one can determine from there 
the distance to the cluster, and thus the Hubble rate Hq. 



The advantages of this method are that it can be apphed to large distances and it is based on 
clear physical principles. The main systematics come from possible dumpiness of the gas (which would 
reduce Hq), projection effects (if the clusters are prolate, Hq could be larger), the assumption of hy- 
drostatic equilibrium of the X-ray gas, details of models for the gas and electron densities, and possible 



contaminations from point sources. Present measurements give the value [26] 



Ho = 60± 10 (Icr statistical) ± 20% (systematic) kms^^Mpc"^ , (64) 

compatible with other determinations. A great advantage of this completely new and independent method 
is that nowadays more and more clusters are observed in the X-ray, and soon we will have high-resolution 
2D maps of the SZ decrement from several balloon flights, as well as from future microwave background 
satellites, together with precise X-ray maps and spectra from the Chandra X-ray observatory recently 
launched by NASA, as well as from the European X-ray satellite XMM launched a few months ago by 
ESA, which will deliver orders of magnitude better resolution than the existing Einstein X-ray satellite. 

3.1.3 Cepheid variability 

Cepheids are low-mass variable stars with a period-luminosity relation based on the helium ionization 
cycles inside the star, as it contracts and expands. This time variability can be measured, and the star's 
absolute luminosity determined from the calibrated relationship. From the observed flux one can then 
deduce the luminosity distance, see Eq. (|6T|), and thus the Hubble rate Hq. The Hubble Space Telescope 
(HST) was launched by NASA in 1990 (and repaired in 1993) with the specific project of calibrating the 
extragalactic distance scale and thus determining the Hubble rate with 10% accuracy. The most recent 



results from HST are the following [ ]28| ] 



ffo = 71 ± 4 (random) it 7 (systematic) kms ^Mpc . (65) 

The main source of systematic error is the distance to the Large Magellanic Cloud, which provides the 
fiducial comparison for Cepheids in more distant galaxies. Other systematic uncertainties that affect the 
value of Hq are the internal extinction correction method used, a possible metallicity dependence of the 
Cepheid period-luminosity relation and cluster population incompleteness bias, for a set of 21 galaxies 
within 25 Mpc, and 23 clusters within z ^ 0.03. 

With better telescopes coming up soon, like the Very Large Telescope (VLT) interferometer of the 
European Southern Observatory (ESO) in the Chilean Atacama desert, with 4 synchronized telescopes 
by the year 2005, and the Next Generation Space Telescope (NGST) proposed by NASA for 2008, it is 
expected that much better resolution and therefore accuracy can be obtained for the determination of Hq. 

3.2 The matter content Q.m 

In the 1920s Hubble realized that the so called nebulae were actually distant galaxies very similar to our 
own. Soon afterwards, in 1933, Zwicky found dynamical evidence that there is possibly ten to a hundred 
times more mass in the Coma cluster than contributed by the luminous matter in galaxies p9|]. However, 
it was not until the 1970s that the existence of dark matter began to be taken more seriously. At that time 
there was evidence that rotation curves of galaxies did not fall off with radius and that the dynamical 
mass was increasing with scale from that of individual galaxies up to clusters of galaxies. Since then, 
new possible extra sources to the matter content of the universe have been accumulating: 

r^M = ^B, lum (stars in galaxies) (66) 

+ ^^B.dark (MACHOs?) (67) 

+ ^CDM (weakly interacting : axion, neutralino?) (68) 

+ f^HDM (massive neutrinos?) (69) 



The empirical route to the determination of 0,m is nowadays one of the most diversified of all 
cosmological parameters. The matter content of the universe can be deduced from the mass-to-light ratio 
of various objects in the universe; from the rotation curves of galaxies; from microlensing and the direct 
search of Massive Compact Halo Objects (MACHOs); from the cluster velocity dispersion with the use 
of the Virial theorem; from the baryon fraction in the X-ray gas of clusters; from weak gravitational 
lensing; from the observed matter distribution of the universe via its power spectrum; from the cluster 
abundance and its evolution; from direct detection of massive neutrinos at SuperKamiokande; from direct 
detection of Weakly Interacting Massive Particles (WIMPs) at DAMA and UKDMC, and finally from 
microwave background anisotropics. I will review here just a few of them. 



3.2.1 Luminous matter 

The most straight forward method of estimating Qm is to measure the luminosity of stars in galaxies and 
then estimate the mass-to-light ratio, defined as the mass per luminosity density observed from an object, 
T = M. / C. This ratio is usually expressed in solar units, M-q/Cq, so that for the sun Tq = 1. The 
luminosity of stars depends very sensitively on their mass and stage of evolution. The mass-to-light ratio 
of stars in the solar neighbourhood is of order T w 3. For globular clusters and spiral galaxies we can 
determine their mass and luminosity independently and this gives T ^ few. For our galaxy, 

£gai = (1.0 ± 0.3) X 10^ h Lq Mpc~^ and Tgai = 6 ± 3 . (70) 

The contribution of galaxies to the luminosity density of the universe (in the visible-V spectral band, 
centered at ~ 5500 A) is Q 

^v- = (l-7±0.6) X lO^/^L0Mpc-^ (71) 
which can be translated into a mass density by multiplying by the observed T in that band, 

nuh = {6.1 ±2.2) X 10"'^ Tv . (72) 
All the luminous matter in the universe, from galaxies, clusters of galaxies, etc., account for T w 10, 



and thus [31] 



0.002 < Jlium h < 0.006 . (73) 

As a consequence, the luminous matter alone is far from the critical density. Moreover, comparing with 
the amount of baryons from Big Bang nucleosynthesis (^), we conclude that f^ium ^ ^^B^ so there must 
be a large fraction of baryons that are dark, perhaps in the form of very dim stars. 



3.2.2 Rotation curves of spiral galaxies 

The flat rotation curves of spiral galaxies provide the most direct evidence for the existence of large 
amounts of dark matter. Spiral galaxies consist of a central bulge and a very thin disk, stabilized against 
gravitational collapse by angular momentum conservation, and surrounded by an approximately spher- 
ical halo of dark matter. One can measure the orbital velocities of objects orbiting around the disk as a 
function of radius from the Doppler shifts of their spectral lines. The rotation curve of the Andromeda 
galaxy was first measured by Babcock in 1938, from the stars in the disk. Later it became possible to 
measure galactic rotation curves far out into the disk, and a trend was found [32]. The orbital velocity 
rose linearly from the center outward until it reached a typical value of 200 km/s, and then remained flat 
out to the largest measured radii. This was completely unexpected since the observed surface luminosity 
of the disk falls off exponentially with radius, I{r) = /q exp(— r/r/)), see Ref. [^. Therefore, one 
would expect that most of the galactic mass is concentrated within a few disk lengths rn, such that the 
rotation velocity is determined as in a Keplerian orbit, Vj-ot = {GM/rY^'^ oc r~^/^. No such behaviour 
is observed. In fact, the most convincing observations come from radio emission (from the 21 cm line) of 



neutral hydrogen in the disk, which has been measured to much larger galactic radii than optical tracers. 
A typical case is that of the spiral galaxy NGC 6503, where = 1.73 kpc, while the furthest measured 
hydrogen line is at r = 22.22 kpc, about 13 disk lengths away. The measured rotation curve is shown in 
Fig. 13 together with the relative components associated with the disk, the halo and the gas. 




Fig. 13: The rotation curve of the spiral galaxy NGC 6503, determined by radio observations of hydrogen gas in the disk p>3j\. 
The dashed line shows the rotation curve expected from the disk material alone, the dot-dashed line is from the dark matter halo 
alone. 

Nowadays, thousands of galactic rotation curves are known, and all suggest the existence of about 
ten times more mass in the halos of spiral galaxies than in the stars of the disk. Recent numerical simula- 



tions of galaxy formation in a CDM cosmology [ ]34| ] suggest that galaxies probably formed by the infall 



of material in an overdense region of the universe that had decoupled from the overall expansion. The 
dark matter is supposed to undergo violent relaxation and create a virialized system, i.e. in hydrostatic 
equilibrium. This picture has led to a simple model of dark-matter halos as isothermal spheres, with 
density profile p{r) = pc/ (r^ + r^), where Vc is a core radius and pc = v'^/A-kG, with Voc, equal to the 
plateau value of the flat rotation curve. This model is consistent with the universal rotation curve seen in 



Fig. [13[ At large radii the dark matter distribution leads to a flat rotation curve. Adding up all the matter 



in galactic halos up to maximum radii, one finds Thaio ^ 30 h, and therefore 

f^haio > 0.03 - 0.05 . (74) 

Of course, it would be extraordinary if we could confirm, through direct detection, the existence of 
dark matter in our own galaxy. For that purpose, one should measure its rotation curve, which is much 
more difficult because of obscuration by dust in the disk, as well as problems with the determination of 
reliable galactocentric distances for the tracers. Nevertheless, the rotation curve of the Milky Way has 
been measured and conforms to the usual picture, with a plateau value of the rotation velocity of 220 
km/s, see Ref. [^]. For dark matter searches, the crucial quantity is the dark matter density in the solar 
neighbourhood, which turns out to be (within a factor of two uncertainty depending on the halo model) 
Pdm = 0.3 GeV/cm^. We will come back to direct searched of dark matter in a later subsection. 

3.2.3 Microlensing 

The existence of large amounts of dark matter in the universe, and in our own galaxy in particular, is now 
established beyond any reasonable doubt, but its nature remains a mystery. We have seen that baryons 
cannot account for the whole matter content of the universe; however, since the contribution of the halo 



([74|) is comparable in magnitude to the baryon fraction of the universe (38), one may ask whether the 
galactic halo could be made of purely baryonic material in some non-luminous form, and if so, how one 
should search for it. In other words, are MACHOs the non-luminous baryons filling the gap between 
^lum and r^B? If not, what are they? 

Let us start a systematic search for possibilities. They cannot be normal stars since they would 
be luminous; neither hot gas since it would shine; nor cold gas since it would absorb light and reemit 
in the infrared. Could they be bumt-out stellar remnants? This seems implausible since they would 
arise from a population of normal stars of which there is no trace in the halo. Neutron stars or black 
holes would typically arise from Supernova explosions and thus eject heavy elements into the galaxy, 
while the overproduction of helium in the halo is strongly constrained. They could be white dwarfs, i.e. 
stars not massive enough to reach supernova phase. Despite some recent arguments, a halo composed 
by white dwarfs is not rigorously excluded. Are they stars too small to shine? Perhaps M-dwarfs, stars 
with a mass M < 0.1 Mq which are intrinsically dim; however, very long exposure images of the 
Hubble Space Telescope restrict the possible M-dwarf contribution to the galaxy to be below 6%. The 
most plausible alternative is a halo composed of brown dwarfs with mass M < 0.08 Mq, which never 
ignite hydrogen and thus shine only from the residual energy due to gravitational contraction.]^ In fact, 
the extrapolation of the stellar mass function to small masses predicts a large number of brown dwarfs 
within normal stellar populations. A final possibility is primordial black holes (PBH), which could have 



been created in the early universe from early phase transitions []36|], even before baryons were formed, 
and thus may be classified as non-baryonic. They could make a large contribution towards the total ^m, 
and still be compatible with Big Bang nucleosynthesis. 
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Fig. 14: Geometry of the light deflection by a pointlike mass which gives two images of a source viewed by an observer. From 
Ref. 0]. 



Whatever the arguments for or against baryonic objects as galactic dark matter, nothing would be 
more convincing than a direct detection of the various candidates, or their exclusion, in a direct search 
experiment. Fortunately, in 1986 Paczynski proposed a method for detecting faint stars in the halo of our 
galaxy [|^]. The idea is based on the well known effect that a point-like mass deflector placed between 
an observer and a light source creates two different images, as shown in Fig. 14. When the source is 
exactly aligned with the deflector of mass Mo, the image would be an annulus, an Einstein ring, with 
radius 

r| = AGMod, where d = J^l^ (75) 

di + d2 

is the reduced distance to the source, see Fig. 14 . If the two images cannot be separated because their 
angular distance a is below the resolving power of the observer's telescope, the only effect will be an 

"a sometimes discussed alternative, planet-size Jupiters, can be classified as low-mass brown dwarfs. 
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Fig. 15: The apparent lightcurve of a source if a pointlike MACHO passes through the line of sight with a transverse velocity v 
and an impact parameter b. The amplification factor A is shown in logarithmic scale to give the usual astronomical magnitude 
of an object. From Ref. 



apparent brightening of tiie star, an effect known as gravitational microlensing. The, amplification factor 

A = — , , where n = — , (76) 

with r the distance from the line of sight to the deflector. Imagine an observer on Earth watching a distant 
star in the Large Magellanic Cloud (LMC), 50 kpc away. If the galactic halo is filled with MACHOs, 
one of them will occasionally pass near the line of sight and thus cause the image of the background 
star to brighten. If the MACHO moves with velocity v transverse to the line of sight, and if its impact 
parameter, i.e. the minimal distance to the line of sight, is h, then one expects an apparent lightcurve 
as shown in Fig. 15 for different values of h/r-^,. The natural time unit is At = te/u, and the origin 
corresponds to the time of closest approach to the line of sight. 

The probability for a target star to be lensed is independent of the mass of the dark matter ob- 
ject [39, 1^]. For stars in the LMC one finds a probability, i.e. an optical depth for microlensing of 
the galactic halo, of approximately r ~ 10~^. Thus, if one looks simultaneously at several millions of 
stars in the LMC during extended periods of time, one has a good chance of seeing at least a few of 
them brightened by a dark halo object. In order to be sure one has seen a microlensing event one has to 



monitor a large sample of stars long enough to identify the characteristic light curve shown in Fig. [15 



The unequivocal signatures of such an event are the following: it must be a) unique (non-repetitive in 
time); b) time-symmetric; and c) achromatic (because of general covariance). These signatures allow 
one to discriminate against variable stars which constitute the background. The typical duration of the 
light curve is the time it takes a MACHO to cross an Einstein radius. At = r-^/v. If the deflector mass 
is 1 Mq, the average microlensing time will be 3 months, for 10~^ Mq it is 9 days, for 10~^ Mq it is 
I day, and for 10~® Mq it is 2 hours. A characteristic event, of duration 34 days, is shown in Fig. [l^. 

The first microlensing events towards the LMC were reported by the MACHO and EROS collab- 
orations in 1993 [40, 41]. Nowadays, there are 12 candidates towards the LMC, 2 towards the SMC, 
around 40 towards the bulge of our own galaxy, and about 2 towards Andromeda, seen by AGAPE [42], 
with a slightly different technique based on pixel brightening rather than individual stars. Thus, mi- 
crolensing is a well established technique with a rather robust future. In particular, it has allowed the 
MACHO and EROS collaboration to draw exclusion plots for various mass ranges in terms of their max- 
imum allowed halo fraction, see Fig. 17. The MACHO Collaboration conclude in their 5-year analysis. 
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Fig. 16: The best candidate (LMC-1) for microlensing from the MACHO Collaboration in the direction of the Large Mage- 
llanic Cloud. A recent reanalysis of this event suggested an amplification factor ^max ~ 7.20 ± 0.09, with achromaticity 
Arad/Auuc = 1.00 ± 0.05, and a duration of i = 34.8 ± 0.2. From Ref. Q. 



see Ref. [38], that the spatial distribution of events is consistent with an extended lens distribution such 
as Milky Way or LMC halo, consisting partially of compact objects. A maximum likelihood analysis 
gives a MACHO halo fraction of 20% for a typical halo model with a 95% confidence interval of 8% to 
50%. A 100% MACHO halo is ruled out at 95% c.l. for all except their most extreme halo model. The 
most likely MACHO mass is between 0.15 Mq and 0.9 Mq, depending on the halo model. The lower 
mass is characteristic of white dwarfs, but a galactic halo composed primarily of white dwarfs is barely 
compatible with a range of observational constraints. On the other hand, if one wanted to attribute the 
observed events to brown dwarfs, one needs to appeal to a very non-standard density and/or velocity dis- 
tribution of these objects. It is still unclear what sort of objects the microlensing experiments are seeing 
towards the LMC and where the lenses are. Nevertheless, the field is expanding, with several new exper- 
iments already underway, to search for clear signals of parallax, or binary systems, where the degeneracy 
between mass and distance can be resolved. For a discussion of those new results, see Ref. [37]. 



3.2.4 Virial theorem and large scale motion 

Clusters of galaxies are the largest gravitationally bound systems in the universe (superclusters are not 
yet in equilibrium). We know today several thousand clusters; they have typical radii of 1 — 5 Mpc and 
typical masses of 2— 9 x 10^'* Mq. Zwicky noted in 1933 that these systems appear to have large amounts 
of dark matter [p9|]. He used the virial theorem (for a gravitationally bound system in equilibrium), 
2(£'kin) = — (i?grav)> whcrc (E'kin) = \rn{v'^) is the average kinetic energy of one of the bound objects 
(galaxies) of mass m and (-Egrav) = —m{GM/r) is the average gravitational potential energy caused 
by the attraction of the other galaxies. Measuring the velocity dispersion (w^) from the Doppler shifts of 
the spectral lines and estimating the geometrical size of the system gives an estimate of its total mass M. 
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Fig. 17: Likelihood contours for MACHO mass m (in units of solar mass) and halo fraction / for a typical size halo. The plus 
sign shows the maximum likelihood estimate and the contours enclose regions of 68%, 90%, 95% and 99% probability. The 
panels are labeled according to different sets of selection criteria (A or B), and whether or not an LMC halo with MACHO 
fraction / is included. From Ref. |B8||. 



As Zwicky noted, this virial mass of clusters far exceeds their luminous mass, typically leading to a 
mass-to-light ratio Truster = 200 ± 70. Assuming that the average cluster T is representative of the 
entire universe one finds for the cosmic matter density [ 44 ] 



= 0.24 ± 0.05 (Icr statistical) ± 0.09 (systematic) . (77) 

On scales larger than clusters the motion of galaxies is dominated by the overall cosmic expansion. 
Nevertheless, galaxies exhibit peculiar velocities with respect to the global cosmic flow. For example, 
our Local Group of galaxies is moving with a speed of 627 it 22 km/s relative to the cosmic microwave 
background reference frame, towards the Great Attractor. 

In the context of the standard gravitational instability theory of structure formation, the peculiar 
motions of galaxies are attributed to the action of gravity during the universe evolution, caused by the 
matter density inhomogeneities which give rise to the formation of structure. The observed large-scale 
velocity fields, together with the observed galaxy distributions, can then be translated into a measure for 
the mass-to-light ratio required to explain the large-scale flows. An example of the reconstruction of the 



matter density field in our cosmological vicinity from the observed velocity field is shown in Fig. [18 



The cosmic matter density inferred from such analyses is |43, 45] 

> 0.3 95% c.l. (78) 
Related methods that are more model-dependent give even larger estimates. 

3.2.5 Baryon fraction in clusters 

Since large clusters of galaxies form through gravitational collapse, they scoop up mass over a large 
volume of space, and therefore the ratio of baryons over the total matter in the cluster should be rep- 
resentative of the entire universe, at least within a 20% systematic error. Since the 1960s, when X-ray 
telescopes became available, it is known that galaxy clusters are the most powerful X-ray sources in the 
sky p6|]. The emission extends over the whole cluster and reveals the existence of a hot plasma with 
temperature T ~ 10^ — 10® K, where X-rays are produced by electron bremsstrahlung. Assuming the 

'^Recent observations indicate that T is independent of scale up to supercluster scales ~ 100 h^^ Mpc. 





Fig. 18: The velocity and density fluctuation fields in the Supergalactic Plane as recovered by the POTENT method from the 
Mark III velocities of about 3,000 galaxies with 12 h^^ smoothing. The vectors are projections of the 3D velocity field in 
the frame of the CMB. Coordinates are in units of 10 h"^ Mpc. The marked structures are the Local Group (LG), the "Great 
Attractor" (GA), the Coma cluster "Great Wall" (GW), the Perseus-Pisces (PP) region and the "Southern Wall" (SW). From 
Ref. 0]. 

gas to be in hydrostatic equilibrium and applying the virial theorem one can estimate the total mass in 
the cluster, giving general agreement (within a factor of 2) with the virial mass estimates. From these 
estimates one can calculate the baryon fraction of clusters 

/b/i3/2 ^ Q Q3 _ Q Qg ^ i^B _ Q ^ for = 0.65 , (79) 

which together with ([73|) indicates that clusters contain far more baryonic matter in the form of hot gas 
than in the form of stars in galaxies. Assuming this fraction to be representative of the entire universe, 
and using the Big Bang nucleosynthesis value of Qb = 0.05 ± 0.01, for h = 0.65, we find 

Qm = 0.3 ±0.1 (statistical) ± 20% (systematic) . (80) 

This value is consistent with previous determinations of Qm- If some baryons are ejected from the cluster 
during gravitational collapse, or some are actually bound in nonluminous objects like planets, then the 
actual value of Qm is smaller than this estimate. 



3.2.6 Weak gravitational lensing 

Since the mid 1980s, deep surveys with powerful telescopes have observed huge arc-like features in 



galaxy clusters, see for instance Fig. |19[ The spectroscopic analysis showed that the cluster and the 
giant arcs were at very different redshifts. The usual interpretation is that the arc is the image of a 
distant background galaxy which is in the same line of sight as the cluster so that it appears distorted 
and magnified by the gravitational lens effect: the giant arcs are essentially partial Einstein rings. From 
a systematic study of the cluster mass distribution one can reconstruct the shear field responsible for the 
gravitational distortion, see Ref. |47]. 




Fig. 19: The most famous image of weak gravitational lensing around the Abell 2218 cluster, made by the Hubble Space 
Telescope. From Ref. [p^. 

This analysis shows that there are large amounts of dark matter in the clusters, in rough agreement 
with the virial mass estimates, although the lensing masses tend to be systematically larger. At present, 
the estimates indicate JIm = 0.2 — 0.3 on scales ^ 6 Mpc, while JIm = 0.4 for the Corona Borealis 
supercluster, on scales of order 20 Mpc. 



3.2. 7 Structure formation and the matter power spectrum 

One the most important constraints on the amount of matter in the universe comes from the present distri- 
bution of galaxies. As we mentioned in the Section 2.3, gravitational instability increases the primordial 
density contrast, seen at the last scattering surface as temperature anisotropics, into the present density 
field responsible for the large and the small scale structure. 

Since the primordial spectrum is very approximately represented by a scale-invariant Gaussian 
random field, the best way to present the results of structure formation is by working with the 2-point 
correlation function in Fourier space (the equivalent to the Green's function in QFT), the so-called power 
spectrum. If the reprocessed spectrum of inhomogeneities remains Gaussian, the power spectrum is all 
we need to describe the galaxy distribution. Non-Gaussian effects are expected to arise from the non- 
linear gravitational collapse of structure, and may be important at small scales [|T^]. 

The power spectrum measures the degree of inhomogeneity in the mass distribution on different 
scales. It depends upon a few basic ingredientes: a) the primordial spectrum of inhomogeneities, whether 
they are Gaussian or non-Gaussian, whether adiabatic (perturbations in the energy density) or isocur- 
vature (perturbations in the entropy density), whether the primordial spectrum has tilt (deviations from 
scale-invariance), etc.; b) the recent creation of inhomogeneities, whether cosmic strings or some other 
topological defect from an early phase transition are responsible for the formation of structure today; and 
c) the cosmic evolution of the inhomogeneity, whether the universe has been dominated by cold or hot 
dark matter or by a cosmological constant since the beginning of structure formation, and also depending 
on the rate of expansion of the universe. 



The working tools used for the comparison between the observed power spectrum and the pre- 
dicted one are very precise N-body numerical simulations and theoretical models that predict the shape 
but not the amplitude of the present power spectrum. Even though a large amount of work has gone 
into those analyses, we still have large uncertainties about the nature and amount of matter necessary for 
structure formation. A model that has become a working paradigm is a flat cold dark matter model with 
a cosmological constant and = 0.3 — 0.4. This model will soon be confronted with very precise 
measurements from SDSS, 2dF, and several other large redshift catalogs, that are already taking data, see 
Section 4.5. 

The observational constraints on the power spectrum have a huge lever arm of measurements 
at very different scales, mainly from the observed cluster abundance, on 10 Mpc scales, to the CMB 
fluctuations, on 1000 Mpc scales, which determines the normalization of the spectrum. At present, deep 
redshift surveys are probing scales between 100 and 1000 Mpc, which should begin to see the turnover 
corresponding to the peak of the power spectrum at k^c^, see Figs. ^ and ^. The standard CDM model with 
r^M = 1> normalized to the CMB fluctuations on large scales, is inconsistent with the cluster abundance. 
The power spectra of both a flat model with a cosmological constant or an open universe with JIm = 0.3 
(defined as ACDM and OCDM, respectively) can be normalized so that they agree with both the CMB 
and cluster observations. In the near future, galaxy survey observations will greatly improve the power 
spectrum constraints and will allow a measurement of I^m from the shape of the spectrum. At present, 
these measurements suggest a low value of Om, but with large uncertainties. 

3.2.8 Cluster abundance and evolution 

Rich clusters are the most recently formed gravitationally bound systems in the universe. Their number 
density as a function of time (or redshift) helps determine the amount of dark matter. The observed 
present (z ~ 0) cluster abundance provides a strong constraint on the normalization of the power spec- 
trum of density perturbations on cluster scales. Both ACDM and OCDM are consistent with the observed 
cluster abundance at z ~ 0, see Fig. while Standard CDM (Einstein-De Sitter model, with f^M = 1)> 
when normalized at COBE scales, produces too many clusters at all redshifts. 
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Fig. 20: The evolution of the cluster abundance as a function of redshift, compared with observations from massive clusters. 
The four models are normalized to COBE. From Ref. l|49||. 



The evolution of the cluster abundance with redshift breaks the degeneracy among the models 
at z ~ 0. The low-mass models (Open and A-CDM) predict a relatively small change in the number 
density of rich clusters as a function of redshift because, due to the low density, hardly any structure 



growth occurs since z ~ 1. The high-mass models (Tilted and Standard CDM) predict that structure has 
grown steadily and rich clusters only formed recently: the number density of rich clusters at z 1 is 
predicted to be exponentially smaller than today. The observation of a single massive cluster is enough to 
rule out the Qm = 1 model. In fact, three clusters have been seen, suggesting a low density universe []50|], 

nu = 0.25 toil (IfT statistical) ± 20% (systematic) . (81) 

But one should be cautious. There is the caveat that for this constraint it is assumed that the initial 
spectrum of density perturbations is Gaussian, as predicted in the simplest models of inflation, but that 
has not yet been confirmed observationally on cluster scales. 
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Fig. 21: The observed cosmic matter components as functions of the Hubble expansion parameter. The luminous matter 
component is given by Eq. (|73|); the galactic halo component is the horizontal band, Eq. ([t^), crossing the baryonic component 
from BBN, Eq. (||); and the dynamical mass component from large scale structure analysis is given by Eq. (|80|). Note that in 
the range Ho = 70 ± 7 km/s/Mpc, there are three dark matter problems, see the text. 



3.2.9 Summary of the matter content 

We can summarize the present situation with Fig. for J^m as a function of Hq. There are four bands, 
the luminous matter f^ium; the baryon content JIb. from BBN; the galactic halo component Ohaio> and 
the dynamical mass from clusters, Q.y{- From this figure it is clear that there are in fact three dark matter 
problems: The first one is where are 90% of the baryons. Between the fraction predicted by BBN and 
that seen in stars and diffuse gas there is a huge fraction which is in the form of dark baryons. They could 
be in small clumps of hydrogen that have not started thermonuclear reactions and perhaps constitute the 
dark matter of spiral galaxies' halos. Note that although $7b and ilhaio coincide at Hq ~ 70 km/s/Mpc, 
this could be just a coincidence. The second problem is what constitutes 90% of matter, from BBN 
baryons to the mass inferred from cluster dynamics. This is the standard dark matter problem and could 
be solved by direct detection of a weakly interacting massive particle in the laboratory. And finally, 
since we know from observations of the CMB, see Section 4.4, that the universe is flat, what constitutes 
around 60% of the energy density, from dynamical mass to critical density, JIq = 1? One possibility 
could be that the universe is dominated by a diffuse vacuum energy, i.e. a cosmological constant, which 
only affects the very large scales. Alternatively, the theory of gravity (general relativity) may need to 
be modified on large scales, e.g. due to quantum gravity effects. The need to introduce an effective 
cosmological constant on large scales is nowadays the only reason why gravity may need to be modified 
at the quantum level. Since we still do not have a quantum theory of gravity, such a proposal is still very 
speculative, and most of the approaches simply consider the inclusion of a cosmological constant as a 
phenomenological parameter. 



3.2.10 Massive neutrinos 



One of the 'usual suspects' when addressing the problem of dark matter are neutrinos. They are the only 
candidates known to exist. If neutrinos have a mass, could they constitute the missing matter? We know 
from the Big Bang theory, see Section 2.2.2, that there is a cosmic neutrino background at a temperature 
of approximately 2K. This allows one to compute the present number density in the form of neutrinos, 
which turns out to be, for massless neutrinos, ny{Ty) = ^ n^{Tj) = 112 cm^'^, per species of neutrino. 



If neutrinos have mass, as recent experiments seem to suggest, see Fig. |22|, the cosmic energy density in 
massive neutrinos would be = J2 n^mi, = Y^rij J2 f^v^ ^^id therefore its contribution today. 



94 eV 



(82) 



The discussion in the previous Sections suggest that ^ 0-4, and thus, for any of the three families 
of neutrinos, rriy < 40 eV. Note that this limit improves by six orders of magnitude the present bound 



on the tau-neutrino mass ||5 1|]. Supposing that the missing mass in non-baryonic cold dark matter arises 
from a single particle dark matter (PDM) component, its contribution to the critical density is bounded 



by 0.05 < QpDMh^ < 0.4, see Fig. |21 
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Fig. 22: The neutrino parameter space, mixing angle against Am^, including the results from the different solar and atmospheric 
neutrino oscillation experiments. Note the threshold of cosmologically important masses, cosmologically detectable neutrinos 
(by CMB and LSS observations), and cosmologically excluded range of masses. From Ref. |p2|]. 



I will now go through the various logical arguments that exclude neutrinos as the dominant compo- 
nent of the missing dark matter in the universe. Is it possible that neutrinos with a mass 4 eV < m,^ < 40 



eV be the non-baryonic PDM component? For instance, could massive neutrinos constitute the dark 
matter halos of galaxies? For neutrinos to be gravitationally bound to galaxies it is necessary that their 
velocity be less that the escape velocity Vesc, and thus their maximum momentum is Pmax = it^u Vqsc- 
How many neutrinos can be packed in the halo of a galaxy? Due to the Pauli exclusion principle, 
the maximum number density is given by that of a completely degenerate Fermi gas with momen- 
tum pf = Pmax, i-c. njnax = Pmax/^TT^- Therefore, the maximum local density in dark matter 
neutrinos is pmax = "-max^ijy = f^sf./37r^, which must be greater than the typical halo density 
Phaio = 0.3 GeVcm"^. For a typical spiral galaxy, this constraint, known as the Tremaine-Gunn limit, 
gives > 40 eV, see Ref. [53]. However, this mass, even for a single species, say the tau-neutrino, 
gives a value for Q^h'^ = 0.5, which is far too high for structure formation. Neutrinos of such a low mass 
would constitute a relativistic hot dark matter component, which would wash-out structure below the su- 
percluster scale, against evidence from present observations, see Fig. Furthermore, applying the same 
phase-space argument to the neutrinos as dark matter in the halo of dwarf galaxies gives rrij^ > 100 eV, 
beyond closure density (p2|). We must conclude that the simple idea that light neutrinos could constitute 
the particle dark matter on all scales is ruled out. They could, however, still play a role as a sub-dominant 
hot dark matter component in a flat CDM model. In that case, a neutrino mass of order 1 eV is not 
cosmological excluded, see Fig. 22 . 

Another possibility is that neutrinos have a large mass, of order a few GeV. In that case, their num- 
ber density at decoupling, see Section 2.2.2, is suppressed by a Boltzmann factor, ~ exp(— mj^/Tdoc). 
For masses m,^ > Tdec — 0.8 MeV, the present energy density has to be computed as a solution 
of the corresponding Boltzmann equation. Apart from a logarithmic correction, one finds Qi^h"^ ~ 
0.1(10 GeV/m,!,)^ for Majorana neutrinos and slightly smaller for Dirac neutrinos. In either case, neu- 
trinos could be the dark matter only if their mass was a few GeV. Laboratory limits for Vj- of around 18 
MeV [51], and much more stringent ones for and Ve, exclude the known light neutrinos. However, 
there is always the possibility of a fourth unknown heavy and stable (perhaps sterile) neutrino. If it 
couples to the Z boson and has a mass below 45 GeV for Dirac neutrinos (39.5 GeV for Majorana neu- 
trinos), then it is ruled out by measurements at LEP of the invisible width of the Z. There are two logical 
alternatives, either it is a sterile neutrino (it does not couple to the Z), or it does couple but has a larger 
mass. In the case of a Majorana neutrino (its own antiparticle), their abundance, for this mass range, 
is too small for being cosmologically relevant, ilj^/i^ < 0.005. If it were a Dirac neutrino there could 
be a lepton asymmetry, which may provide a higher abundance (similar to the case of baryogenesis). 
However, neutrinos scatter on nucleons via the weak axial-vector current (spin-dependent) interaction. 
For the small momentum transfers imparted by galactic WIMPs, such collisions are essentially coherent 
over an entire nucleus, leading to an enhancement of the effective cross section. The relatively large 
detection rate in this case allowes one to exclude fourth-generation Dirac neutrinos for the galactic dark 
matter f^^. Anyway, it would be very implausible to have such a massive neutrino today, since it would 
have to be stable, with a life-time greater than the age of the universe, and there is no theoretical reason 
to expect a massive sterile neutrino that does not oscillate into the other neutrinos. 

Of course, the definitive test to the possible contribution of neutrinos to the overall density of 
the universe would be to measure directly their mass in laboratory experiments. There are at present 
two types of experiments: neutrino oscillation experiments, which measure only differences in squared 
masses, and direct mass-searches experiments, like the tritium /3-spectrum and the neutrinoless double-/? 
decay experiments, which measure directly the mass of the electron neutrino and give a bound m^^ ^ 2 
eV. Neutrinos with such a mass could very well constitute the HDM component of the universe, r^HDM ~ 
0.15. The oscillation experiments give a variety of possibilities for Am^ = 0.3 — 3 eV^ from LSND 
(not yet confirmed), to the atmospheric neutrino oscillations from SuperKamiokande (Amf, ~ 3 x 
10^'^ eV^) and the solar neutrino oscillations (A ~ 10~^ eV^). Only the first two possibilities would 
be cosmologically relevant, see Fig. 22. 



'^For a review of Neutrinos, see Bilenky's contribution to these Proceedings [p5[ 



3.2.11 Weakly Interacting Massive Particles 

Unless we drastically change the theory of gravity on large scales, baryons cannot make up the bulk 
of the dark matter. Massive neutrinos are the only alternative among the known particles, but they are 
essentially ruled out as a universal dark matter candidate, even if they may play a subdominant role as 
a hot dark matter component. There remains the mystery of what is the physical nature of the dominant 
cold dark matter component. 

Something like a heavy stable neutrino, a generic Weakly Interacting Massive Particle (WIMP), 
could be a reasonable candidate because its present abundance could fall within the expected range. 



3 X 10^27 pj^^Sg 

(o'ann'Wrel) 



(83) 



Here v^d is the relative velocity of the two incoming dark matter particles and the brackets (...) denote 
a thermal average at the freeze-out temperature, Tf ~ mpDM/20, when the dark matter particles go out 
of equilibrium with radiation. The value of (o"ann^^rci) needed for fipDM ~ 1 is remarkably close to 
what one would expect for a WIMP with a mass jtipdm — 100 GeV, (cann^rei) 



a 



■in mpDM 



3 X 10 cm s . We still do not know whether this is just a coincidence or an important hint on the 
nature of dark matter. 
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Fig. 23: The maximum likelihood region from the annual-modulation signal consistent with a neutralino of mass — 59 1]^4 



GeV and a proton cross section of ^cTp = 7.0 l"'2 x 10 pb, see the text. The scatter plot represents the theoretical predictions 
of a generic MSSM. From Ref. 



There are a few theoretical candidates for WIMPs, like the neutralino, coming from supersymme- 
tric extensions of the standard model of particle physics,[3 but at present there is no empirical evidence 
that such extensions are indeed realized in nature. In fact, the non-observation of supersymmetric par- 
ticles at current accelerators places stringent limits on the neutralino mass and interaction cross sec- 



tion [57]. 



If WIMPs constitute the dominant component of the halo of our galaxy, it is expected that some 
may cross the Earth at a reasonable rate to be detected. The direct experimental search for them rely 



For a review of Supersymmetry (SUSY), see Carena's contribution to these Proceedings. 



on elastic WIMP collisions with the nuclei of a suitable target. Dark matter WIMPs move at a typical 
galactic virial velocity of around 200 — 300 km/s, depending on the model. If their mass is in the 
range 10 — 100 GeV, the recoil energy of the nuclei in the elastic collision would be of order 10 keV. 
Therefore, one should be able to identify such energy depositions in a macroscopic sample of the target. 
There are at present three different methods: First, one could search for scintillation light in Nal crystals 
or in liquid xenon; second, search for an ionization signal in a semiconductor, typically a very pure 
germanium crystal; and third, use a cryogenic detector at 10 mK and search for a measurable temperature 
increase of the sample. The main problem with such a type of experiment is the low expected signal rate, 
with a typical number below 1 event/kg/day. To reduce natural radioactive contamination one must 
use extremely pure substances, and to reduce the background caused by cosmic rays requires that these 
experiments be located deeply underground. 
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Fig. 24: The DAMA experiment sees an annual variation, of order 7%, in the WIMP flux due to the Earth's motion around the 
Sun. The model independent residual rate in the lowest (2 — 6 keV) cumulative energy interval (in counts per day/kg/keV) is 
shown as a function of time since 1 January of the first year of data taking. The expected behaviour of a WIMP signal is a 
cosine function with a minimum (maximum) roughly at the dashed (dotted) vertical lines. From Ref. nSq]. 



The best limits on WIMP scattering cross sections come from some germanium experiments [58], 
as well as from the Nal scintillation detectors of the UK dark matter collaboration (UKDMC) in the 



Boulby salt mine in England []59|], and the DAMA experiment in the Gran Sasso laboratory in Italy 
Current experiments already touch the parameter space expected from supersymmetric particles, see 
Fig. 23, and therefore there is a chance that they actually discover the nature of the missing dark matter. 
The problem, of course, is to attribute a tentative signal unambiguously to galactic WIMPs rather than to 
some unidentified radioactive background. 

One specific signature is the annual modulation which arises as the Earth moves around the Sun.[^ 
Therefore, the net speed of the Earth relative to the galactic dark matter halo varies, causing a modulation 
of the expected counting rate. The DAMA/Nal experiment has actually reported such a modulation sig- 
nal, see Fig. from the combined analysis of their 4-year data [|^, which provides a confidence level of 
99.6% for a neutralino mass of = 52 Ig*^ GeV and a proton cross section of ^cjp = 7.2 [l^Q g x 10~^ 
pb, where ^ = p^(/0.3 GeVcm^'^ is the local neutralino energy density in units of the galactic halo 
density. There has been no confirmation yet of this result from other dark matter search groups, but 
hopefully in the near future we will have much better sensitivity at low masses from the Cryogenic Rare 
Event Search with Superconducting Thermometers (CRESST) experiment at Gran Sasso as well as at 
weaker cross sections from the CDMS experiment at Stanford and the Soudan mine, see Fig. 25. The 
CRESST experiment [50] uses sapphire crystals as targets and a new method to simultaneously measure 

'^The time scale of the Sun's orbit around the center of the galaxy is too large to be relevant in the analysis. 




Fig. 25: Exclusion range for the spin-independent WIMP scattering cross section per nucleon from the Nal experiments and 
the Ge detectors. Also shown is the range of expected counting rates for neutralinos in the MSSM. The search goals for the 
upcoming large-scale cryogenic detectors CRESST and CDMS are also shown. From Ref. n23]. 



the phonons and the scintillating light from particle interactions inside the crystal, which allows excel- 
lent background discrimination. Very recently there has been the interesting proposal of a completely 
new method based on a Superheated Droplet Detector (SDD), which claims to have already a similar 
sensitivity as the more standard methods described above, see Ref. [jsl]]. 

There exist other indirect methods to search for galactic WIMPs [^]. Such particles could self- 
annihilate at a certain rate in the galactic halo, producing a potentially detectable background of high 
energy photons or antiprotons. The absence of such a background in both gamma ray satellites and the 
Alpha Matter Spectrometer |63] imposes bounds on their density in the halo. Alternatively, WIMPs 
traversing the solar system may interact with the matter that makes up the Earth or the Sun so that a 
small fraction of them will lose energy and be trapped in their cores, building up over the age of the 
universe. Their annihilation in the core would thus produce high energy neutrinos from the center of the 
Earth or from the Sun which are detectable by neutrino telescopes. In fact. Super Kamiokande already 
covers a large part of SUSY parameter space. In other words, neutrino telescopes are already competitive 
with direct search experiments. In particular, the AMANDA experiment at the South Pole [Q], which 
is expected to have 10^ Cherenkov detectors 2.3 km deep in very clear ice, over a volume ~ 1 km'^, is 
competitive with the best direct searches proposed. The advantages of AMANDA are also directional, 
since the arrays of Cherenkov detectors will allow one to reconstruct the neutrino trajectory and thus its 
source, whether it comes from the Earth or the Sun. 



3.3 The cosmological constant JIa 

A cosmological constant is a term in the Einstein equations, see Eq. ([T]), that corresponds to the energy 
density of the vacuum of quantum field theories, A = SirGpy, see Ref. []65|]. These theories predict 
a value of order ~ Mp ~ 5 x 10^^ g/cm^, which is about 123 orders of magnitude larger than 
the critical density ([l^. Such a discrepancy is one of the biggest problems of theoretical physics []6^. 
It has always been assumed that quantum gravity effects, via some as yet unknown symmetry, would 
exactly cancel the cosmological constant, but this remains a downright speculation. Moreover, one of the 



difficulties with a non-zero value for A is that it appears coincidental that we are now living at a special 
epoch when the cosmological constant starts to dominate the dynamics of the universe, and that it will do 
so forever after, see Section 2.1.2 and Eq. (20). Nevertheless, ever since Einstein introduced it in 1917, 
this ethereal constant has been invoked several times in history to explain a number of apparent crises, 
always to disappear- under further scrutiny [21 1. 
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Fig. 26: The Type la supernovae observed nearby show a relationship between their absolute luminosity and the timescale of 
their light curve: the brighter supernovae are slower and the fainter ones are faster. A simple linear relation between the absolute 
magnitude and a "stretch factor" multiplying the light curve timescale fits the data quite well. From Ref. ||68|]. 



In spite of the theoretical prejudice towards A = 0, there are new observational arguments for a 
non-zero value. The most compelling ones are recent evidence that we live in a flat universe, from obser- 
vations of CMB anisotropics, together with strong indications of a low mass density universe {Q-m < 1)> 
from the large scale distribution of galaxies, clusters and voids, that indicate that some kind of dark 
energy must make up the rest of the energy density up to critical, i.e. = 1 — JIm- In addition, the 
discrepancy between the ages of globular clusters and the expansion age of the universe may be cleanly 
resolved with A 7^ 0. Finally, there is growing evidence for an accelerating universe from observations 
of distant supernovae. I will now discuss the different arguments one by one. 

The only known way to reconcile a low mass density with a flat universe is if an additional "dark" 
energy dominates the universe today. It would have to resist gravitational collapse, otherwise it would 
have been detected already as part of the energy in the halos of galaxies. However, if most of the energy 
of the universe resists gravitational collapse, it is impossible for structure in the universe to grow. This 
dilemma can be resolved if the hypothetical dark energy was negligible in the past and only recently 
became the dominant component. According to general relativity, this requires that the dark energy have 
negative pressure, since the ratio of dark energy to matter density goes like a{t)~^P/ ^. This argument [57] 
would rule out almost all of the usual suspects, such as cold dark matter, neutrinos, radiation, and kinetic 



energy, since they all have zero or positive pressure. Thus, we expect something like a cosmological 
constant, with negative pressure, p ^ —p,to account for the missing energy. 

This negative pressure would help accelerate the universe and reconcile the expansion age of the 
universe with the ages of stars in globular clusters, see Fig. |n|, where toHo is shown as a function 
of r2M> in a flat universe, Qa = 1 — ^M, and an open one, JIa = 0. For the present age of the 
universe of to = 13 ± 1 Gyr, and the measured rate of expansion, Hq = 70 ± 7 km/s/Mpc, one finds 



to^^o = 0.93 ± 0.12 (adding errors in quadrature), which corresponds to Qm = 0.05 



+0.24 
-0.10 



for an open 

universe, see Fig. |11), marginally consistent with observations of large scale structure. On the other hand, 
for a flat universe with a cosmological constant, to^^o = 0.93 ± 0.12 corresponds to = 0.34 11q'i2, 
which is perfectly compatible with recent observations. These suggest that we probably live in a flat 
universe that is accelerating, dominated today by a vacuum energy density. 
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Fig. 27: Hubble diagram for the high redshift supernovae found by the SN Cosmology Project. From Ref. [|68|]. A similar 
diagram is found by the High Redshift Supernova Project [^. Both groups conclude that distant supernovae are fainter than 
expected, and this could be due to an accelerating universe. 



This conclusions have been supported by growingly robust observational evidence from distant 
supernovae. In their quest for the cosmological parameters, astronomers look for distant astrophysical 
objects that can serve as standard candles to determine the distance to the object from their observed 
apparent luminosity. A candidate that has recently been exploited with great success is a certain type of 
supernova explosions at large redshifts, called SN of type la. These are white dwarf stars at the end of 
their life cycle that accrete matter from a companion until they become unstable and violently explode 
in a natural thermonuclear explosion that out-shines their progenitor galaxy. The intensity of the distant 
flash varies in time, it takes about three weeks to reach its maximum brightness and then it declines over 
a period of months. Although the maximum luminosity varies from one supernova to another, depending 
on their original mass, their environment, etc., there is a pattern: brighter explosions last longer than 
fainter ones. By studying the characteristic light curves, see Fig. of a reasonably large statistical 
sample, cosmologists from two competing groups, the Supernova Cosmology Project [^] and the High- 



redshift Supernova Project [69], are confident that they can use this type of supernova as a standard 
candle. Since the light coming from some of these rare explosions has travelled for a large fraction of 
the size of the universe, one expects to be able to infer from their distribution the spatial curvature and 
the rate of expansion of the universe. 

One of the surprises revealed by these observations is that high redshift type la supernovae appear 
fainter than expected for either an open (Om < 1) or a flat (Qm = 1) universe, see Fig. 27. In fact, the 
universe appears to be accelerating instead of decelerating, as was expected from the general attraction 




Fig. 28: The left figure shows the best-fit confidence regions (68% - 99% c.l.) in the (S7m, ^Ia) plane, for the high redshift 
supemovae results. The systematic uncertainty is not shown, and would shift the ellipses vertically. Present observations 
disfavour the Eisntein-de Sitter model (circle) by several standard deviations. The upper-left shaded region represents "bouncing 
universe" cosmologies with no Big Bang in the past. The lower-right shaded region corresponds to a universe that is younger 
than the oldest heavy elements, for any value of h > 0.5. The right figure shows the isochrones of constant Hoto, the age of 
the universe in units of the Hubble time, H^^, with the best-fit 68% and 90% confidence regions in the (flm, ^a) plane. From 
Ref. ||6§]. 

of matter, see Eq. (|22|); something seems to be acting as a repulsive force on very large scales. The 
most natural explanation for this is the presence of a cosmological constant, a diffuse vacuum energy 
that permeates all space and, as explained above, gives the universe an acceleration that tends to separate 
gravitationally bound systems from each other. The best-fit results from the Supernova Cosmology 
Project give a linear combination 0.8Om— O-Gf^A = — 0.2ib0.1 (Icj), and, for a flat universe (J7m+^^a = 
1), the best-fit values for the combined analysis of both groups [l68|, [69|], are 

= 0.28 to^ol {la statistical) (identified systematics) , (84) 

nl^^ = 0.72 to^ol {la statistical) +[];|]^ (identified systematics) . (85) 

However, one may think that it is still premature to conclude that the universe is indeed accel- 
erating, because of possibly large systematic errors inherent to most cosmological measurements, and 
in particular to observations of supemovae at large redshifts. There has been attempts to find crucial 
systematic effects like evolution, chemical composition dependence, reddening by dust, etc. in the su- 
pemovae observations that would invalidate the claims, but none of them are now considered as a serious 
threat. Perhaps the most critical one today seems to be sampling effects, since the luminosities of the 
high-redshift supemovae (z ~ 0.5 — 1.0) are all measured relative to the same set of local supemovae 
{z < 0.3). Hence, absolute calibrations, completeness levels, and any other systematic effects related 
to both data sets are critical. For instance, the intense efforts to search for high-redshift objects have 
led to the peculiar situation where the nearby sample, which is used for calibration, is now smaller than 
the distant one. Further searches, already underway, for increasing the nearby supemovae sample will 
provide an important check. 

Moreover, there are bounds on a cosmological constant that come from the statistics of gravita- 
tional lensing, with two different methods. Gravitational lensing can be due to various accumulations of 
matter along the line of sight to the distant light sources. The first method uses the abundance of multiply 
imaged sources like quasars, lensed by intervening galaxies. The probability of finding a lensed image 



is directly proportional to the number of galaxies (lenses) along the path and thus to the distance to the 
source. This distance, for fixed Hq, increases dramatically for a large value of the cosmological constant: 
the age of the universe and the distance to the galaxy become large for 7^ because the universe has 
been expanding for a longer time; therefore, more lenses are predicted for Q,/^ > 0. Using this method, 
an upper limit of 

Qa < 0.75 (95% c.l.) (86) 



has recently been obtained [JTOp, marginally consistent with the supernovae results, but there are caveats 
to this powerful method due to uncertainties in the number density and lensing cross section of the lensing 
galaxies as well as the distant quasars. A second method is lensing by massive clusters of galaxies, which 
produces widely separated lensed images of quasars and distorted images of background galaxies. The 
observed statistics, when compared with numerical simulations, rule out the Qm = 1 models and set an 



upper bound on the cosmological constant, JIa < 0.7, see Ref. [47]. However, this limit is very sensitive 



to the resolution of the numerical simulations, which are currently improving. 
3.4 The spatial curvature ^Ik 

As we will discuss in detail in Section 4.4, observations of the two-point correlation function of tem- 
perature anisotropics in the microwave background provide a crucial test for the spatial curvature of the 
universe. From those observations one can tell whether the photons that left the last scattering surface, 
at redshift z = 1100, have travelled in straight lines, like in a flat universe, or in curved paths, like in 
an open one. Very recent observations made by the balloon experiment BOOMERANG suggest that the 
universe is indeed spatially flat (Qk = 0) with about 10% accuracy ^n\\, 

17o = jIm + ^^a = 1-0±0.1 (95% c.l.) (87) 

These measuremnts are bound to be improved in the near future, by both balloon experiments and by the 
Microwave Anisotropy Probe (MAP) satellite, to be launched by NASA at the end of year 2000 



Furthermore, with the launch in 2007 of Planck satellite p3\\ we will be able to determine Qq with 1% 
accuracy. 

3.5 The age of the universe to 

The universe must be older than the oldest objects it contains. Those are believed to be the stars in the 
oldest clusters in the Milky Way, globular clusters. The most reliable ages come from the application 
of theoretical models of stellar evolution to observations of old stars in globular clusters. For about 30 



years, the ages of globular clusters have remained reasonable stable, at about 15 Gyr [74]. However, 



recently these ages have been revised downward [75]. 



During the 1980s and 1990s, the globular cluster age estimates have improved as both new obser- 
vations have been made with CCDs, and since refinements to stellar evolution models, including opaci- 
ties, consideration of mixing, and different chemical abundances have been incorporated [}76|]. From the 
theory side, uncertainties in globular cluster ages come from uncertainties in convection models, opac- 
ities, and nuclear reaction rates. From the observational side, uncertainties arise due to corrections for 
dust and chemical composition. However, the dominant source of systematic errors in the globular clus- 
ter age is the uncertainty in the cluster distances. Fortunately, the Hipparcos satellite recently provided 
geometric parallax measurements for many nearby old stars with low metallicity, typical of glubular clus- 
ters, thus allowing for a new calibration of the ages of stars in globular clusters, leading to a downward 
revision to 10 — 13 Gyr [[7^]. Moreover, there were very few stars in the Hipparcos catalog with both 
small parallax erros and low metal abundance. Hence, an increase in the sample size could be critical in 
reducing the statatistical uncertaintites for the calibration of the globular cluster ages. There are already 
proposed two new parallax satellites, NASA's Space Interferometry Mission (SIM) and ESA's mission, 
called GAIA, that will give 2 or 3 orders of magnitude more accurate parallaxes than Hipparcos, down 



to fainter magnitude limits, for several orders of magnitude more stars. Until larger samples are avail- 
able, however, distance errors are likely to be the largest source of systematic uncertainty to the globular 
cluster age [21]. 
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Fig. 29: The recent estimates of the age of the universe and that of the oldest objects in our galaxy. The last three points 
correspond to the combined analysis of 8 different measurements, for h — 0.64, 0.68 and 7.2, which indicates a relatively weak 
dependence on h. The age of the Sun is accurately known and is included for reference. Error bars indicate la limits. The 
averages of the ages of the Galactic Halo and Disk are shaded in gray. Note that there isn't a single age estimate more than 
2a away from the average. The result to > ^gai is logically inevitable, but the standard EdS model does not satisfy this unless 
h < 0.55. FromRef. 

The supernovae groups can also determine the age of the universe from their high redshift ob- 



servations. Figure 28 shows that the confidence regions in the (r^Mj^^A) plane are almost parallel to 
the contours of constant age. For any value of the Hubble constant less than Hq = 70 km/s/Mpc, the 
implied age of the universe is greater than 13 Gyr, allowing enough time for the oldest stars in globular 
clusters to evolve ^J^. Integrating over JIm and Oa, the best fit value of the age in Hubble-time units is 
Hoto = 0.93 ± 0.06 or equivalently to = 14.1 ± 1.0 (0.65 h'^) Gyr \E^]. The age would be somewhat 



larger in a fiat universe: i^o^o'^* = 0-96 Vqqj or, equivalently. 



+0.09 



.flat 
'^O 



14.4 t\i 



(0.65 h^') Gyr 



(88) 



Furthermore, a combination of 8 independent recent measurements: CMB anisotropics, type la 
SNe, cluster mass-to-light ratios, cluster abundance evolution, cluster baryon fraction, deuterium-to- 
hidrogen ratios in quasar spectra, double-lobed radio sources and the Hubble constant, can be used to 
determine the present age of the universe |77]. The result is shown in Fig. compared to other recent 
determinations. The best fit value for the age of the universe is, according to this analysis, to = 13.4ibl.6 
Gyr, about a billion years younger than other recent estimates 



We can summarize this Section by showing the region in parameter space where we stand nowa- 



days, thanks to the recent cosmological observations. We have plotted that region in Fig. 30. One could 
also superimpose the contour lines corresponding to equal to^^o lines, as a cross check. It is extraordi- 
nary that only in the last few months we have been able to reduce the concordance region to where it 
stands today, where all the different observations seem to converge. There are still many uncertainties, 
mainly systematic; however, those are quickly decreasing and becoming predominantly statistical. In the 
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Fig. 30: The concordance region. The sum Hm + gives the total cosmic energy content and determines the geometry of 
space-time, whether spatially flat, open or closed. Their difference, Q,m /2 — Q,a , characterizes the relative strength of expansion 
and gravity, and determines how the expansion rate changes with time, whether accelerating or decelerating. Furthermore, a 
balance between the two densities determines the fate of the universe, whether it will expand forever or recoUapse. These three 
effects have been probed by recent observations, from large scale structure (cluster data), temperature anisotropics (microwave 
background data) and the universe expansion (supernova data). Surprisingly enough, at present all observations seem to lie 
within a naroow region in parameter space. The Einstein-de Sitter model is no longer the preferred one. The best model today 
is a flat model with a third of the energy density in the form of non-relativistic matter and two thirds in the form of vacuum 
energy or a cosmological constant. From Ref. j?^]. 

near future, with precise observations of the anisotropics in the microwave background temperature and 
polarization, to be discussed in Section 4.4, we will be able to reduce those uncertainties to the level of 
one percent. This is the reason why cosmologists are so excited and why it is claimed that we live in the 
Golden Age of Cosmology. 

4 THE INFLATIONARY PARADIGM 

The hot Big Bang theory is nowadays a very robust edifice, with many independent observational checks: 
the expansion of the universe; the abundance of light elements; the cosmic microwave background; a 
predicted age of the universe compatible with the age of the oldest objects in it, and the formation of 
structure via gravitational collapse of initially small inhomogeneities. Today, these observations are 
confirmed to within a few percent accuracy, and have helped establish the hot Big Bang as the preferred 
model of the universe. All the physics involved in the above observations is routinely tested in the 
laboratory (atomic and nuclear physics experiments) or in the solar system (general relativity). 

However, this theory leaves a range of crucial questions unanswered, most of which are initial 
conditions' problems. There is the reasonable assumption that these cosmological problems will be 
solved or explained by new physical principles at high energies, in the early universe. This assumption 
leads to the natural conclusion that accurate observations of the present state of the universe may shed 
light onto processes and physical laws at energies above those reachable by particle accelerators, present 
or future. We will see that this is a very optimistic approach indeed, and that there are many unresolved 
issues related to those problems. However, there might be in the near future reasons to be optimistic. 



4.1 Shortcomings of Big Bang Cosmology 

The Big Bang theory could not explain the origin of matter and structure in the universe; that is, the 
origin of the matter-antimatter asymmetry, without which the universe today would be filled by a uniform 
radiation continuosly expanding and cooling, with no traces of matter, and thus without the possibility 
to form gravitationally bound systems like galaxies, stars and planets that could sustain life. Moreover, 
the standard Big Bang theory assumes, but cannot explain, the origin of the extraordinary smoothness 
and flatness of the universe on the very large scales seen by the microwave background probes and the 
largest galaxy catalogs. It cannot explain the origin of the primordial density perturbations that gave rise 
to cosmic structures like galaxies, clusters and superclusters, via gravitational collapse; the quantity and 
nature of the dark matter that we believe holds the universe together; nor the origin of the Big Bang itself. 

A summary [79] of the problems that the Big Bang theory cannot explain is: 

• The global structure of the universe. 

- Why is the universe so close to spatial flatness? 

- Why is matter so homogeneously distributed on large scales? 

• The origin of structure in the universe. 

- How did the primordial spectrum of density perturbations originate? 

• The origin of matter and radiation. 

- Where does all the energy in the universe come from? 

- What is the nature of the dark matter in the universe? 

- How did the matter-antimatter asymmetry arise? 

• The initial singularity. 

- Did the universe have a beginning? 

- What is the global structure of the universe beyond our observable patch? 
Let me discuss one by one the different issues: 



4.1.1 The Flatness Problem 

The Big Bang theory assumes but cannot explain the extraordinary spatial flatness of our local patch of 
the universe. In the general FRW metric @ the parameter K that characterizes spatial curvature is a free 
parameter. There is nothing in the theory that determines this parameter a priori. However, it is directly 
related, via the Friedmann equation to the dynamics, and thus the matter content, of the universe, 

K = ^pa--H-a- = ^pa^(^) • (89) 
We can therefore define a new variable, 

— 1 const. 



Q pa 



2 



(90) 



whose time evolution is given by 
where N = ln(a/aj) characterizes the number of e-folds of universe expansion (dN = Hdt) and 



X = ^ = {l + ?,uj)x, (91) 



where we have used Eq. ( [30[ ) for the time evolution of the total energy, pa?, which only depends on the 
barotropic ratio lo. It is clear from Eq. (^Tj) that the phase-space diagram (x, x') presents an unstable 
critical (saddle) point at x = for cj > —1/3, i.e. for the radiation (u = 1/3) and matter {uj = 0) eras. 
A small perturbation from x = will drive the system towards x = itcxD. Since we know the universe 
went through both the radiation era (because of primordial nucleosynthesis) and the matter era (because 
of structure formation), tiny deviations from i7 = 1 would have grown since then, such that today 

Xo = = Xin (p^) (1 + ^eq) ■ (92) 



eq- 



In order that today's value be in the range 0.1 < Qq < 1-2, or xq w it is required that at, say, 

primordial nucleosynthesis (Tj^g ~ 10^ Teq) its value be 

^7(^^s) = l±lo-l^ (93) 

which represents a tremendous finetuning. Perhaps the universe indeed started with such a peculiar 
initial condition, but it is epistemologically more satisfying if we give a fundamental dynamical reason 
for the universe to have started so close to spatial flatness. These arguments were first used by Robert 
Dicke in the 1960s, much before inflation. He argued that the most natural initial condition for the 
spatial curvature should have been the Planck scale curvature, ^^^R = 6K/ /| , where the Planck length 
is Ip = {fiG/c^Y^'^ = 1.62 X lO^'^^ cm, that is, 60 orders of magnitude smaller than the present size 
of the universe, ao = 1.38 x 10^^ cm. A universe with this immense curvature would have collapsed 
within a Planck time, tp = {hG/c'y/^ = 5.39 x 10"'^^ s, again 60 orders of magnitude smaller than 
the present age of the universe, to = 4.1 x 10^^ s. Therefore, the flatness problem is also related to the 
Age Problem, why is it that the universe is so old and flat when, under ordinary circumstances (based on 
the fundamental scale of gravity) it should have lasted only a Planck time and reached a size of order the 
Planck length? As we will see, inflation gives a dynamical reason to such a peculiar initial condition. 



4.1.2 The Homogeneity Problem 

An expanding universe has particle horizons, that is, spatial regions beyond which causal communica- 
tion cannot occur. The horizon distance can be defined as the maximum distance that light could have 
travelled since the origin of the universe 

ft d+i 

dn{t) = a{t) / —-^H-\t), (94) 
Jo a[t') 

which is proportional to the Hubble scale.[^ For instance, at the beginning of nucleosynthesis the horizon 
distance is a few light-seconds, but grows linearly with time and by the end of nucleosynthesis it is a 
few light-minutes, i.e. a factor 100 larger, while the scale factor has increased only a factor of 10. The 
fact that the causal horizon increases faster, dn ~ t, than the scale factor, a ~ t^/^, implies that at any 
given time the universe contains regions within itself that, according to the Big Bang theory, were never 
in causal contact before. For instance, the number of causally disconnected regions at a given redshift z 
present in our causal volume today, dii{to) = ao, is 

which, for the time of decoupling, is of order A'^cD(^dcc) ~ 10^ ^ 1. 

This phenomenon is particularly acute in the case of the observed microwave background. Infor- 
mation cannot travel faster than the speed of light, so the causal region at the time of photon decoupling 
could not be larger than (iH(tdcc) ~ 3 x 10^ light years across, or about 1° projected in the sky today. So 
why should regions that are separated by more than 1° in the sky today have exactly the same tempera- 
ture, to within 10 ppm, when the photons that come from those two distant regions could not have been 



in causal contact when they were emitted? This constitutes the so-called horizon problem, see Fig. 3 1 



and was first discussed by Robert Dicke in the 1970s as a profound inconsistency of the Big Bang theory. 



4.2 Cosmological Inflation 

In the 1980s, a new paradigm, deeply rooted in fundamental physics, was put forward by Alan H. 



Guth [plU, Andrei D. Linde [ ]82p and others \ p3[ |84[ , |85p, to address these fundamental questions. Ac- 
cording to the inflationary paradigm, the early universe went through a period of exponential expansion. 



For the radiation era, tlie iiorizon distance is equal to tiie Hubble scale. For the matter era it is twice the Hubble scale. 



Our Hubble 




Fig. 3 1 : Perhaps the most acute problem of the Big Bang theory is explaining the extraordinary homogeneity and isotropy of the 
microwave background, see Fig. ^. At the time of decoupling, the volume that gave rise to our present universe contained many 
causally disconnected regions (top figure). Today we observe a blackbody spectrum of photons coming from those regions and 
they appear to have the same temperature, Ti — T2, to one part in 10'''. Why is the universe so homogeneous? This constitutes 
the so-called horizon problem, which is spectacularly solved by inflation. From Ref. ||8C|, [78|]. 



driven by the approximately constant energy density of a scalar field called the inflaton. In modem 
physics, elementary particles are represented by quantum fields, which resemble the familiar electric, 
magnetic and gravitational fields. A field is simply a function of space and time whose quantum oscil- 
lations are interpreted as particles. In our case, the inflaton field has, associated with it, a large potential 
energy density, which drives the exponential expansion during inflation, see Fig. 32. We know from gen- 
eral relativity that the density of matter determines the expansion of the universe, but a constant energy 
density acts in a very peculiar way: as a repulsive force that makes any two points in space separate at 
exponentially large speeds. (This does not violate the laws of causality because there is no information 
carried along in the expansion, it is simply the stretching of space-time.) 

This superluminal expansion is capable of explaining the large scale homogeneity of our observ- 
able universe and, in particular, why the microwave background looks so isotropic: regions separated 
today by more than 1° in the sky were, in fact, in causal contact before inflation, but were stretched to 
cosmological distances by the expansion. Any inhomogeneities present before the tremendous expansion 
would be washed out. This explains why photons from supposedly causally disconneted regions have 



actually the same spectral distribution with the same temperature, see Fig. 31 



Moreover, in the usual Big Bang scenario a flat universe, one in which the gravitational attraction 
of matter is exactly balanced by the cosmic expansion, is unstable under perturbations: a small deviation 
from flatness is amplified and soon produces either an empty universe or a collapsed one. As we dis- 
cussed above, for the universe to be nearly flat today, it must have been extremely flat at nucleosynthesis, 
deviations not exceeding more than one part in 10^^. This extreme fine tuning of initial conditions was 
also solved by the inflationary paradigm, see Fig. 33. Thus inflation is an extremely elegant hypothesis 
that explains how a region much, much greater that our own observable universe could have become 
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Fig. 32: The inflaton field can be represented as a ball rolling down a hill. During inflation, the energy density is approximately 
constant, driving the tremendous expansion of the universe. When the ball starts to oscillate around the bottom of the hill, 
inflation ends and the inflaton energy decays into particles. In certain cases, the coherent oscillations of the inflaton could 
generate a resonant production of particles which soon thermalize, reheating the universe. From Ref . iTa] . 



smooth and flat without recourse to ad hoc initial conditions. Furthermore, inflation dilutes away any 
"unwanted" relic species that could have remained from early universe phase transitions, like monopoles, 
cosmic strings, etc., which are predicted in grand unified theories and whose energy density could be so 
large that the universe would have become unstable, and collapsed, long ago. These relics are diluted by 
the superluminal expansion, which leaves at most one of these particles per causal horizon, making them 
harmless to the subsequent evolution of the universe. 

The only thing we know about this peculiar scalar field, the inflaton, is that it has a mass and 
a self-interaction potential V{(j)) but we ignore everything else, even the scale at which its dynamics 
determines the superluminal expansion. In particular, we still do not know the nature of the inflaton field 
itself, is it some new fiindamental scalar field in the electroweak symmetry breaking sector, or is it just 
some effective description of a more fundamental high energy interaction? Hopefully, in the near future, 
experiments in particle physics might give us a clue to its nature. Inflation had its original inspiration in 
the Higgs field, the scalar field supposed to be responsible for the masses of elementary particles (quarks 
and leptons) and the breaking of the electroweak symmetry. Such a field has not been found yet, and its 
discovery at the future particle colliders would help understand one of the truly fundamental problems in 
physics, the origin of masses. If the experiments discover something completely new and unexpected, it 
would automatically affect the idea of inflation at a fundamental level. 

4. 2. 1 Homogeneous scalar field dynamics 

In this subsection I will describe the theoretical basis for the phenomenon of inflation. Consider a scalar 
field (p, a singlet under any given interaction, with an effective potential V{cl)). The Lagrangian for such 
a field in a curved background is 



£inf = l9'"'d^(t^d,4>-V{cl>), 



(96) 



Fig. 33: The exponential expansion during inflation made the radius of curvature of the universe so large that our observable 
patch of the universe today appears essentialy flat, analogous (in three dimensions) to how the surface of a balloon appears 
flatter and flatter as we inflate it to enormous sizes. This is a crucial prediction of cosmological inflation that will be tested to 
extraordinary accuracy in the next few years. From Ref. [^4[[7^. 



whose evolution equation in a Friedmann-Robeitson- Walker metric @ and for a homogeneous field (f){t) 
is given by 

+ 3F,^ + y' (</.) = 0, (97) 
where H is the rate of expansion, together with the Einstein equations, 

= ^{\^'' + y{<P))^ (98) 



3 V2 

H = -y0', (99) 

where = SttG. The dynamics of inflation can be described as a perfect fluid with a time dependent 
pressure and energy density given by 

P = \4>'' + v{(t>), (100) 

P = \^''-v{4>)- (101) 

The field evolution equation ( |97| ) can then be written as the energy conservation equation, 

p + ?,H{p + p) = . (102) 

If the potential energy density of the scalar field dominates the kinetic energy, V{<f)) ^ 0^, then we see 
that 

p c:^ —p ^ p ~ const. ^ H{(f)) ~ const. , (103) 



which leads to the solution 



a{t) ~ exp(Ht) ^ — > accelerated expansion . (104) 

a 

Using the definition of the number of e-folds, N = ln(a/aj), we see that the scale factor grows expo- 
nentially, a{N) = Ojexp(A^). This solution of the Einstein equations solves immediately the flatness 
problem. Recall that the problem with the radiation and matter eras is that 17 = 1 (x = 0) is an unsta- 
ble critical point in phase-space. However, during inflation, with p ~ —p =^ ~ — 1, we have that 
1 + 3u; > and therefore x = is a stable attractor of the equations of motion, see Eq. (^T|). As a con- 
sequence, what seemed an ad hoc initial condition, becomes a natural prediction of inflation. Suppose 
that during inflation the scale factor increased e-folds, then 

:ro = Xine-2^(^)'(l + Zeq) =^ e-2^10^'^<l ^ iV > 65 , (105) 

-'eq 

where we have assumed that inflation ended at the scale T4nd> and the transfer of the inflaton energy 
density to thermal radiation at reheating occurred almost instantaneously]^ at the temperature Trh ~ 
Knd ~ 10^^ GeV. Note that we can now have initial conditions with a large uncertainty, — 1, and 
still have today xq ^ 1, thanks to the inflationary attractor towards = 1. This can be understood very 
easily by realizing that the three curvature evolves during inflation as 

(3)/? = ^ = , for AT » 1 . (106) 

Therefore, if cosmological inflation lasted over 65 e-folds, as most models predict, then today the uni- 



verse (or at least our local patch) should be exactly flat, see Fig. |33|, a prediction that can be tested with 
great accuracy in the near future and for which already seems to be some evidence from observations of 
the microwave background [|7l|]. 

Furthermore, inflation also solves the homogeneity problem in a spectacular way. First of all, due 
to the superluminal expansion, any inhomogeneity existing prior to inflation will be washed out, 

f— j oc e-2^ 0, foriV»l. (107) 

Moreover, since the scale factor grows exponentially, while the horizon distance remains essentially 
constant, dH{t) ~ H^^ = const., any scale within the horizon during inflation will be stretched by the 
superluminal expansion to enormous distances, in such a way that at photon decoupling all the causally 
disconnected regions that encompass our present horizon actually come from a single region during 
inflation, about 65 e-folds before the end. This is the reason why two points separated more than 1° 
in the sky have the same backbody temperature, as observed by the COBE satellite: they were actually 
in causal contact during inflation. There is at present no other proposal known that could solve the 
homogeneity problem without invoquing an acausal mechanism like inflation. 

Finally, any relic particle species (relativistic or not) existing prior to inflation will be diluted by 
the expansion, 

PM OC ~ e"^^ — > , for » 1 , (108) 

PR OC a-^ ~ e"^^ — > , for > 1 . (109) 

Note that the vacuum energy density p^ remains constant under the expansion, and therefore, very soon 
it is the only energy density remaining to drive the expansion of the universe. 



"There could be a small delay in thermalization, due to the intrinsic inefficiency of reheating, but this does not change 
significantly the required number of e-folds. 



4.2.2 The slow-roll approximation 

In order to simplify the evolution equations during inflation, we will consider the slow-roll approximation 
(SRA). Suppose that, during inflation, the scalar field evolves very slowly down its effective potential, 
then we can define the slow -roll parameters [[8^, 

< 1, (110) 



H_ _ 



6 = --^ < 1. (Ill) 
H(t) 



It is easy to see that the condition 



e<l ^ ->0 (112) 

a 

characterizes inflation: it is all you need for superluminal expansion, i.e. for the horizon distance to grow 
more slowly than the scale factor, in order to solve the homogeneity problem, as well as for the spatial 
curvature to decay faster than usual, in order to solve the flatness problem. 



The number of e-folds during inflation can be written with the help of Eq. (110) as 



N = In^ = / Hdt = ^=^, (113) 



which is an exact expression in terms of e((/)). 



In the limit given by Eqs. (UlQ), the evolution equations (97) and (Bq) become 



h\i-'-) ^ = yy(c/>), (114) 

^ 3H^ = -¥'{(!,). (115) 

Note that this corresponds to a reduction of the dimensionality of phase-space from two to one dimen- 
sions, H{(j),(j)) — > H{(j)). In fact, it is possible to prove a theorem, for single-field inflation, which 
states that the slow-roll approximation is an attractor of the equations of motion, and thus we can al- 
ways evaluate the inflationary trajectory in phase-space within the SRA, therefore reducing the number 
of initial conditions to just one, the initial value of the scalar field. If H{(j)) only depends on 0, then 
= —k'^4>/2 and we can rewrite the slow-roll parameters ( ]11(]| ) as 

(116) 

= r/ - e < 1. (117) 



K? H{(j)) K? V{(p) 2k2 V V{<p) J 



The last expression defines the new slow-roll parameter rj, not to be confused with conformal time (see 
next Section). The number of e-folds can also be rewritten in this approximation as 



a very useful expression for evaluating N for a given effective scalar potential V{(t)). 



4.3 The origin of density perturbations 

If cosmological inflation made the universe so extremely flat and homogeneous, where did the galaxies 
and clusters of galaxies come from? One of the most astonishing predictions of inflation, one that was not 
even expected, is that quantum fluctuations of the inflaton field are stretched by the exponential expansion 
and generate large-scale perturbations in the metric. Inflaton fluctuations are small wave packets of 
energy that, according to general relativity, modify the space-time fabric, creating a whole spectrum of 
curvature perturbations. The use of the word spectrum here is closely related to the case of light waves 
propagating in a medium: a spectrum characterizes the amplitude of each given wavelength. In the 
case of inflation, the inflaton fluctuations induce waves in the space-time metric that can be decomposed 
into different wavelengths, all with approximately the same amplitude, that is, corresponding to a scale- 
invariant spectrum. These patterns of perturbations in the metric are like fingerprints that unequivocally 
characterize a period of inflation. When matter fell in the troughs of these waves, it created density 
perturbations that collapsed gravitationally to form galaxies, clusters and superclusters of galaxies, with 
a spectrum that is also scale invariant. Such a type of spectrum was proposed in the early 1970s (before 
inflation) by Harrison and Zel'dovich [ 17], to explain the distribution of galaxies and clusters of galaxies 
on very large scales in our observable universe. Perhaps the most interesting aspect of structure formation 
is the possibility that the detailed knowledge of what seeded galaxies and clusters of galaxies will allow 
us to test the idea of inflation. 



4.3.1 Gauge invariant perturbation theory 

Until now we have considered only the unperturbed FRW metric described by a scale factor a{t) and a 
homogeneous scalar field 4>{t), 

ds^ = a^{rj)[drj'^ --fijdx'dx^], (119) 
</> = Hv), (120) 

where rj = J dt/a{t) is the conformal time, under which the background equations of motion can be 
written as 

= y Q<^' V aV(0)) , (121) 

n' -n^ = -i^(t>^\ (122) 

0" + 2?^(/.' + oV((/)) = 0, (123) 

where H = aH and (j)' = acf). 

During inflation, the quantum fluctuations of the scalar field will induce metric perturbations which 
will backreact on the scalar field. Let us consider, in linear perturbation theory, the most general line ele- 
ment with both scalar and tensor metric perturbations [[87|],P| together with the scalar field perturbations 

ds^ = a^{r}) [(1 + 2A)drf - 2B\idx'dr] - {(1 + 2n)-iij + 2E\ij + 2hij^dx'dx^^ , (124) 

= ct){r,) + 6^{v,x'). (125) 

The indices {i,j} label the three-dimensional spatial coordinates with metric jij, and the \i denotes 
covariant derivative with respect to that metric. The gauge invariant tensor perturbation hij corresponds 
to a transverse traceless gravitational wave, V^hij = h} = 0. The four scalar perturbations {A, B,Tl, E) 
are gauge dependent functions of (r/, x*). Under a general coordinate (gauge) transformation [|^, Q 

f) = r? + e°(r?,:E*), (126) 
ii = +-f'^i\Ari,x'), {111) 



'*Note that inflation cannot generate, to linear order, a vector perturbation. 



with arbitrary functions {^^,(,), the scalar and tensor perturbations transform, to linear order, as 



'jz = n-nf, 

hi 



B = B + f-C', 
E = E-C, 



(128) 
(129) 
(130) 



where a prime denotes derivative with respect to conformal time. It is possible to construct, however. 



(131) 
(132) 



two gauge-invariant gravitational potentials [K7, 88], 



^ = A+{B- E')' + nifi - E') , 

^ = n + n{B - E'), 



which are related through the perturbed Einstein equations. 



-dp, 



(133) 
(134) 



where 5 pis the gauge-invariant density perturbation, and the latter expression is nothing but the Poisson 
equation for the gravitational potential, written in relativistic form. 

During inflation, the energy density is given in terms of a scalar field, and thus the gauge-invariant 
equations for the perturbations on comoving hypersurfaces (constant energy density hypersurfaces) are 



(135) 
(136) 
(137) 

6(1)" + 2n6(l)' -V^6(l) = 4(/.'$'-2aV(</.)$-aV'(0)50. (138) 

This system of equations seem too difficult to solve at first sight. However, there is a gauge 
invariant combination of variables that allows one to find exact solutions. Let us define [88] 

u = a5(j) + z<^, (139) 

z^a^. (140) 
n 

Under this redefinition, the above equations simplify enormously to just three independent equations, 

(141) 



z 



u" - V^n u = , 

z 

V'^ = !f-li(zn'-z'u), 
2 



-zu . 



(142) 
(143) 



From Equation ( 141 ) we can find a solution u{z), which substituted into (143) can be integrated to give 
^{z), and together with u{z) allow us to obtain 6(p{z). 



4.3.2 Quantum Field Theory in curved space-time 

Until now we have treated the perturbations as classical, but we should in fact consider the perturbations 
$ and 5(f) as quantum fields. Note that the perturbed action for the scalar mode u can be written as 



6S 



1 



d^xdr] \{u'f - {Vuf + —u" 



(144) 



In order to quantize the field u in the curved background defined by the metric (1 19), we can write the 
operator 



u(r/,x) 



(145) 



(2^)3/2 

where the creation and annihilation operators satisfy the commutation relation of bosonic fields, and the 
scalar field's Fock space is defined through the vacuum condition, 

[ak,4'] = '^'(k-k'), (146) 
OklO) = 0. (147) 

Note that we are not assuming that the inflaton is a fundamental scalar field, but that is can be written as 
a quantum field with its commutation relations (as much as a pion can be described as a quantum field). 

The equations of motion for each mode Uk (??) are decoupled in linear perturbation theory. 



(148) 



The ratio z" / z acts like a time-dependent potential for this Schrodinger like equation. In order to find 
exact solutions to the mode equation, we will use the slow-roll parameters (110), see Ref. |86] 



1 



1 



H(t)' 



K? Z^ 

= 1 + e 



z 

TTz 



(149) 
(150) 



In terms of these parameters, the conformal time and the effective potential for the Uk mode can be 
written as 



-1 



+ 



eda 



^=n\l + e-5){2-5)+ n-\e' - S')] . 
z 

Note that the slow -roll parameters, (149) and ( |150D , can be taken as constant]^ to order e^, 

e' = 2neie-5) = 0{e^), 



In that case, for constant parameters, we can write 



z 
z 



1/2 1\ 



H4>' 

_ -1 1 

where v 



l+e-5 1 



(151) 
(152) 

(153) 
(154) 

(155) 
(156) 



"For instance, there are models of inflation, like power-law inflation, a{t) ~ t*", where e = S — 1/p < 1, that give constant 
slow-roll parameters. 



We are now going to search for approximate solutions of the mode equation (148), where the 
effective potential (152) is of order z" jz ~ TH? in the slow-roll approximation. In quasi-de Sitter there 
is a characteristic scale given by the (event) horizon size or Hubble scale during inflation, . There 
will be modes with physical wavelengths much smaller than this scale, kja S> R, that are well 
within the de Sitter horizon and therefore do not feel the curvature of space-time. On the other hand, 
there will be modes with physical wavelengths much greater than the Hubble scale, k/a <§:i H. In these 
two asymptotic regimes, the solutions can be written as 



1 



Uk = Ciz 



-ikrj 



k <^aH . 



(157) 
(158) 



In the limit k » aH the modes behave like ordinary quantum modes in Minkowsky space-time, ap- 
propriately normalized, while in the opposite limit, u/z becomes constant on superhorizon scales. For 



approximately constant slow -roll parameters one can find exact solutions to (148), with the effective 



potential given by (156), that interpolate between the two asymptotic solutions. 



(159) 



where Hy \z) is the Hankel function of the first kind [p9|], and v is given by (156) in terms of the 
slow-roll parameters. In the limit ki] —>■ 0, the solution becomes 



\Uk\ 



_ r>) 
2k r(|) 



V2k ^aHJ 



C{v) 

We can now compute $ and 



r>) 
r(i) 



-ki]Y 



{l-eY-2 ~ 1 for e,(5< 1, 



(160) 



(161) 



from the super-Hubble-scale mode solution (158), for k <^ aH 



Substituting into Eq. (143), we find 



Ci 



n 

7^2 



a drj 5" . 



n 



(162) 
(163) 



The term proportional to Ci corresponds to the growing solution, while that proportional to C2 corre- 
sponds to the decaying solution, which can soon be ignored. These quantities are gauge invariant but 
evolve with time outside the horizon, during inflation, and before entering again the horizon during the 
radiation or matter eras. We would like to write an expression for a gauge invariant quantity that is 
also constant for superhorizon modes. Fortunately, in the case of adiabatic perturbations, there is such a 
quantity: 

1 u 
C = $ + — + = -, (164) 



which is constant, see Eq. ( |158| ), for k <^ aH. In fact, this quantity C is identical, for superhorizon 
modes, to the gauge invariant curvature metric perturbation TZc on comoving (constant energy density) 



hypersurfaces, see Ref. [ ]87| , |90p . 



1 



eW2 



(165) 



Using Eq. (142) we can write the evolution equation for ( = ^ as ^' ~ 7h which confirms that 
( is constant for (adiabatic^) superhorizon modes, k ^ aH. Therefore, we can evaluate the Newtonian 



This conservation fails for entropy or isocurvature perturbations, see Ref. 



potential when the perturbation reenters the horizon during radiation/matter eras in terms of the 
curvature perturbation TZf^ when it left the Hubble scale during inflation, 



1 



n 



a dr] IZk 



3 + 3a; 
5 + 3a; 



radiation era, 
matter era . 



(166) 



Let us now compute the tensor or gravitational wave metric perturbations generated during infla- 
tion. The perturbed action for the tensor mode can be written as 



5S 



i2 , 



d X dr] 



2k2 



ij) 



with the tensor field hij considered as a quantum field, 



f d^k r 
J E [hk{v)eij{KX)ak,xe'^-'' + h.c. 



(167) 



(168) 



(2vr)3/2 ^^^^^ 

where Cjj (k, A) are the two polarization tensors, satisfying symmetric, transverse and traceless conditions 



k e,; 



0. 



0, 



e,,(-k, A) = e:,(k, A) , ^ e^(k, A)e'^ (k, A) = 4 . 



(169) 
(170) 



while the creation and annihilation operators satisfy the usual commutation relation of bosonic fields, 
Eq. (146). We can now redefine our gauge invariant tensor amphtude as 



V2t 



hkiv) , 



(171) 



which satisfies the following evolution equation, decoupled for each mode Vk{r]) in linear perturbation 
theory, 

v'^+(k'--)v, = 0. 



(172) 

The ratio a" /a acts like a time-dependent potential for this Schrodinger like equation, analogous to the 
term z" / z for the scalar metric perturbation. For constant slow-roll parameters, the potential becomes 



a V 2 



1 1 



We can solve equation ( |172[ ) in the two asymptotic regimes, 

1 



Vk 



Vk 



-ikrj 



'2k 
Ca 



k:$> aH , 
k <^aH . 



(173) 
(174) 

(175) 
(176) 



In the limit k » aH the modes behave like ordinary quantum modes in Minkowsky space-time, ap- 
propriately normalized, while in the opposite limit, the metric perturbation hk becomes constant on 
superhorizon scales. For constant slow -roll parameters one can find exact solutions to (172), with effec- 
tive potential given by (173), that interpolate between the two asymptotic solutions. These are identical 
to Eq. (159) except for the substitution u fi. In the limit krj —>■ 0, the solution becomes 

V2k ^aHJ 



\Vk\ 



(177) 



Since the mode becomes constant on superhorizon scales, we can evaluate the tensor metric pertur- 
bation when it reentered during the radiation or matter era directly in terms of its value during inflation. 



4.3.3 Power spectrum of scalar and tensor metric perturbations 

Not only do we expect to measure the amplitude of the metric perturbations generated during inflation 
and responsible for the anisotropics in the CMB and density fluctuations in LSS, but we should also be 
able to measure its power spectrum, or two-point correlation function in Fourier space. Let us consider 
first the scalar metric perturbations TZ^, which enter the horizon at a = k/ H. Its correlator is given 
by [|8§] 

m*MQ) = ^ _ k') ^ ^gl (2^)3 53(k _ k') , (178) 

^ /,N Wk\'^ / H\2 ^ k \3-2u ,n / k \n-l 

f«W = ^^ = ^y . (179) 



where we have used TZ^ = Ck = ^ and Eq. ( |160| ). This last equation determines the power spectrum in 
terms of its amplitude at horizon-crossing, As, and a tilt, 

n-l^^^^^^ = 3-2. = 2(^):.2,-6e, (180) 
dink VI — e/ 



see Eqs. ( |116[ ), ( \IIT\ ). Note from this equation that it is possible, in principle, to obtain from inflation a 
scalar tilt which is either positive (n > 1) or negative (n < 1). Furthermore, depending on the particular 
inflationary model [^], we can have significant departures from scale invariance. 

Let us consider now the tensor (gravitational wave) metric perturbation, which enter the horizon 

at a = k/H, 

T.i^\hl,hk>,x\0) = 4 ^ \vk\'6Hk - k') ^ ^ (27r)3 S^k - k') , (181) 

A 



where we have used Eqs. ( |171[ ) and ( |177[ ). Therefore, the power spectrum can be approximated by a 
power-law expression, with amplitude At and tilt 

dlnVJk) „ „ / 2e 

riT -= ~ ■ - — I 



dink 



3-2^i = -(^-^) ~-2e<0, (183) 



which is always negative. In the slow-roll approximation, e ^ 1, the tensor power spectrum is scale 
invariant. 



4.4 The anisotropics of the microwave background 

The metric fluctuations generated during inflation are not only responsible for the density perturbations 
that gave rise to galaxies via gravitational collapse, but one should also expect to see such ripples in 
the metric as temperature anisotropics in the cosmic microwave background, that is, minute deviations 
in the temperature of the blackbody spectrum when we look at different directions in the sky. Such 
anisotropics had been looked for ever since Penzias and Wilson's discovery of the CMB, but had eluded 
all detection, until CORE satellite discovered them in 1992, see Fig. ^ The reason why they took so 
long to be discovered was that they appear as perturbations in temperature of only one part in 10^. Soon 
after COBE, other groups quickly confirmed the detection of temperature anisotropics at around 30 ^K, 
at higher multipole numbers or smaller angular scales. There are at this moment dozens of ground 
and balloon-borne experiments analysing the anisotropics in the microwave background with angular 
resolutions from 10° to a few arc minutes in the sky, see Fig. 34. 



4.4.1 Acoustic oscillations in the plasma 



The physics of the CMB anisotropies is relatively simple ||92|]. The universe just before recombination 
is a very tightly coupled fluid, due to the large electromagnetic Thomson cross section (46). Photons 
scatter off charged particles (protons and electrons), and carry energy, so they feel the gravitational 
potential associated with the perturbations imprinted in the metric during inflation. An overdensity of 
baryons (protons and neutrons) does not collapse under the effect of gravity until it enters the causal 
Hubble radius. The perturbation continues to grow until radiation pressure opposes gravity and sets up 
acoustic oscillations in the plasma, very similar to sound waves. Since overdensities of the same size will 
enter the Hubble radius at the same time, they will oscillate in phase. Moreover, since photons scatter 
off these baryons, the acoustic oscillations occur also in the photon field and induces a pattern of peaks 



in the temperature anisotropies in the sky, at different angular scales, see Fig. [34 



There are three different effects that determine the temperature anisotropies we observe in the 
CMB. First, gravity: photons fall in and escape off gravitational potential wells, characterized by <I> in the 
comoving gauge, and as a consequence their frequency is gravitationally blue- or red-shifted, 5v/v = <I>. 
If the gravitational potential is not constant, the photons will escape from a larger or smaller potential well 
than they fell in, so their frequency is also blue- or red-shifted, a phenomenon known as the Rees-Sciama 
effect. Second, pressure: photons scatter off baryons which fall into gravitational potential wells and the 
two competing forces create acoustic waves of compression and rarefaction. Finally, velocity: baryons 
accelerate as they fall into potential wells. They have minimum velocity at maximum compression and 
rarefaction. That is, their velocity wave is exactly 90° off-phase with the acoustic waves. These waves 
induce a Doppler effect on the frequency of the photons. 

The temperature anisotropy induced by these three effects is therefore given by [ p2| ] 
5T r'^o . 1 (5/9 r • V 

— (r) = $(r,tdec)+2 / ^{r,t)dt + . (184) 

Metric perturbations of different wavelengths enter the horizon at different times. The largest wave- 
lengths, of size comparable to our present horizon, are entering now. There are perturbations with wave- 
lengths comparable to the size of the horizon at the time of last scattering, of projected size about 1° 
in the sky today, which entered precisely at decoupling. And there aie perturbations with wavelengths 
much smaller than the size of the horizon at last scattering, that entered much earlier than decoupling, all 
the way to the time of radiation-matter equality, which have gone through several acoustic oscillations 
before last scattering. All these perturbations of different wavelengths leave their imprint in the CMB 
anisotropies. 

The baryons at the time of decoupling do not feel the gravitational attraction of perturbations with 
wavelength greater than the size of the horizon at last scattering, because of causality. Perturbations with 
exactly that wavelength are undergoing their first contraction, or acoustic compression, at decoupling. 



Those perturbations induce a large peak in the temperature anisotropies power spectrum, see Fig. 34 



Perturbations with wavelengths smaller than these will have gone, after they entered the Hubble scale, 
through a series of acoustic compressions and rarefactions, which can be seen as secondary peaks in 
the power spectrum. Since the surface of last scattering is not a sharp discontinuity, but a region of 
Az ~ 100, see Fig. ^, there will be scales for which photons, travelling from one energy concentration 
to another, will erase the perturbation on that scale, similarly to what neutrinos or HDM do for structure 
on small scales. That is the reason why we don't see all the acoustic oscillations with the same amplitude, 
but in fact they decay exponentialy towards smaller angular scales, an effect known as Silk damping, due 
to photon diffusion [^, p2|]. 



4.4.2 The Sachs-Wolfe effect 



The anisotropies corresponding to large angular scales are only generated via gravitational red-shift 
and density perturbations through the Einstein equations, Sp/p = — 2<I> for adiabatic perturbations; we 
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Fig. 34: There are at present dozens of ground and balloon-borne experiments looking at the microwave background tem- 
perature anisotropies with angular resolutions from 10° to a few arc minutes in the sky, corresponding to multipole numbers 
1 = 2 — 3000. Present observations suggest the existence of a peak in the angular distribution, as predicted by inflation. The 
theoretical curve (thick line) illustrates a particular model which fits the data. From Ref. 



can ignore the Doppler contribution, since the perturbation is non-causal. In that case, the temperature 



anisotropy in the sky today is given by [ |95| ] 



6T 



^$(??Ls)(3(?/0,^, 



+ 2 1 dr'^'{7]o-r)Q{r,e,4>), 

''l-LS 



(185) 



where r]Q is the coordinate distance to the last scattering surface, i.e. the present conformal time, while 
r/LS — determines that comoving hypersurface. The above expression is known as the Sachs-Wolfe 



effect []95|], and contains two parts, the intrinsic and the Integrated Sachs-Wolfe (ISW) effect, due to 
integration along the line of sight of time variations in the gravitational potential. 

In linear perturbation theory, the scalar metric perturbations can be separated into $(r/, x) = 
$(77) Q(x), where <5(x) are the scalar harmonics, eigenfunctions of the Laplacian in three dimensions. 



-k"^ Qkim{f: ^) </')■ These functions have the general form 



(186) 



Qkim{r, 9, (j)) = Iiki{r) Yim{0, 0) , 

where Yim{0, (p) are the usual spherical harmonics []89|]. 

In order to compute the temperature anisotropy associated with the Sachs-Wolfe effect, we have 
to know the evolution of the metric perturbation during the matter era. 



+ 3W + a^A $ - 2K ^> = . 



(187) 



In the case of a flat universe without cosmological constant, the Newtonian potential remains constant 
during the matter era and only the intrinsic SW effect contributes to 6T/T. In case of a non- vanishing 
A, since its contribution is negligible in the past, see Eq. (pi]), most of the photon's trajectory towards us 
is unperturbed, and the only difference with respect to the A = case is an overall factor [66]. We will 
consider here the approximation <I> = const, during the matter era and ignore that factor, see Ref. [97]. 



In a flat universe, the radial part of the eigenfunctions ( |186| ) can be written as [ ]96| ] 



n 



kl[r) 



kji{kr 



(188) 



where ji (z) are the spherical Bessel functions []89|]. The growing mode solution of the metric perturbation 



that left the Hubble scale during inflation contributes to the temperature anisotropics on large scales ( |1 85[ ) 
as 

<5T,„ 1 . . „ 1 



I 



-Y{0,<t>) 



3^(??Ls)(3 



= X! X! 



(189) 



1=2 m=-l 



where we have used the fact that at reentry (at the surface of last scattering) the gauge invariant Newtonian 
potential $ is related to the curvature perturbation TZ at Hubble-crossing during inflation, see Eq. (166); 
and we have expanded 5T /T m spherical harmonics. 

We can now compute the two-point correlation function or angular power spectrum, C{0), of the 
CMB anisotropies on large scales, defined as an expansion in multipole number, 



c{0) 



6T* 6T 
— (n) — (n) 



nn — cos f 
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4tt 



(190) 



1=2 



where Pi{z) are the Legendre polynomials [p9|], and we have averaged over different universe realiza- 



tions. Since the coefficients ai^ are isotropic (to first order), we can compute the Ci = (|a;^ 

dk 



as 
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(S) 



47r 
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k 



rn{k)jf{krio) 



(191) 



where we have used Eqs. (189) and (178). In the case of scalar metric perturbation produced during 



inflation, the scalar power spectrum at reentry is given by Vyzik) = A|(/cryo^ 
approximation, see Eq. (|179|). In that case, one can integrate (|191|) to give 



n-l 



(S) 



2tt 
25 



Ijl + 1) q 

27r 



ni- 

25 



constant , 



n-l i ' 
2 J 



in the power-law 



(192) 



for n = 1 . 



(193) 



This last expression corresponds to what is known as the Sachs-Wolfe plateau, and is the reason why the 
coefficients Ci are always plotted multiplied by /(/ + 1), see Fig. |4[ 

Tensor metric perturbations also contribute with an approximately constant angular power spec- 
trum, /(/ + 1)C;. The Sachs-Wolfe effect for a gauge invariant tensor perturbation is given by [95] 



T 



■no 



dr h'{r]o - r) Qrr{r, 0, < 



where Qrr is the rr-component of the tensor harmonic along the line of sight 
tion h during the matter era satisfies the following evolution equation 

h'l + Sn h'f^ + {k'^ + 2K)hk = 0, 



(194) 

The tensor perturba- 
(195) 



which depends on the wavenumber k, contrary to what happens with the scalar modes, see Eq. (187). For 
a flat {K = 0) universe, the solution to this equation is hk{rj) = h Gk{r]), where h is the constant tensor 
metric perturbation at horizon crossing and Gk{ri) = 3 ji{kr])/kr], normalized so that Gfc(O) = 1 at the 
surface of last scattering. The radial part of the tensor harmonic Qrr in a flat universe can be written 
as ||9^] 



jl - 1)1(1 + + 2) 
nk"^ 



1/2 



(196) 



The tensor angular power spectrum can finally be expressed as 



Cf ) = ^ - !)/(/ + + 2) J^^ ^ V,{k) I, 

j2{xo - x)ji{x) 
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dx '■ 



{xq — 



(197) 
(198) 



where x = krj, and Vg{k) is the primordial tensor spectrum ( |182| ). For a scale invariant spectrum, 
nr = 0, we can integrate (|197|) to give []98|] 



36 



487r" 
"385" 



(199) 



with = (1.1184, 0.8789, . . . , 1.00) for / = 2, 3, ... , 30. Therefore, /(/ + 1) C^^^ also becomes 
constant for large /. Beyond / ~ 30, the Sachs-Wolfe expression is not a good approximation and the 
tensor angular power spectrum decays very quickly at large /, see Fig. 40. 



4.4.3 The consistency relation 

In spite of the success of inflation in predicting a homogeneous and isotropic background on which to 
imprint a scale-invariant spectrum of inhomogeneities, it is difficult to test the idea of inflation. A CMB 
cosmologist before the 1980s would have argued that ad hoc initial conditions could have been at the 
origin of the homogeneity and flatness of the universe on large scales, while a LSS cosmologist would 
have agreed with Harrison and Zel'dovich that the most natural spectrum needed to explain the formation 
of structure was a scale-invariant spectrum. The surprise was that inflation incorporated an understanding 
of both the globally homogeneous and spatially flat background, and the approximately scale-invariant 
spectrum of perturbations in the same formalism. But that could have been a coincidence, and is not 
epistemologically testable. 

What is unique to inflation is the fact that inflation determines not just one but two primordial 
spectra, corresponding to the scalar (density) and tensor (gravitational waves) metric perturbations, from 
a single continuous function, the inflaton potential 1^(0). In the slow-roll approximation, one determines, 
from y ((/>), two continuous functions, Vnik) and Vg{k), that in the power-law approximation reduces 
to two amplitudes. As and At, and two tilts, n and nx- It is clear that there must be a relation between 
the four parameters. Indeed, one can see from Eqs. ( |199 ) and (193) that the ratio of the tensor to scalar 



contribution to the angular power spectrum is proportional to the tensor tilt 

cP 25 / 487r2N 
R^^ = -[l + ^)2e^-2.nT. (200) 

This is a unique prediction of inflation, which could not have been postulated a priori by any cosmol- 
ogist. If we finally observe a tensor spectrum of anisotropics in the CMB, or a stochastic gravitational 



wave background in laser interferometers like LIGO or VIRGO [99], with sufficient accuracy to deter- 
mine their spectral tilt, one might have some chance to test the idea of inflation, via the consistency 
relation ([200|). For the moment, observations of the microwave background anisotropics suggest that 



the Sachs-Wolfe plateau exists, see Fig. ^ but it is still premature to determine the tensor contribution. 
Perhaps in the near future, from the analysis of polarization as well as temperature anisotropics, with the 
CMB satellites MAP and Planck, we might have a chance of determining the validity of the consistency 
relation. 

Assuming that the scalar contribution dominates over the tensor on large scales, i.e. R <^ 1, one 
can actually give a measure of the amplitude of the scalar metric perturbation from the observations of 



the Sachs- Wolfe plateau in the angular power spectrum [|97|], 



iji + 1) q 

27r 



(5)1 1/2 



— = (1.03 ±0.07) X 10"^ 
5 



n 



(201) 
(202) 



1.02 ±0.12. 

These measurements can be used to normalize the primordial spectrum and determine the parameters 



of the model of inflation ||91|]. In the near future these parameters will be determined with much better 
accuracy, as described in Section 4.4.5. 
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Fig. 35: Theoretical predictions for CMB temperature angular power spectra as a function of multipole number / for models 



with primordial adiabatic perturbations. Each graph shows the effect of a variation in one of these parameters. From Ref. [ IOC ]. 



4.4.4 The acoustic peaks 

The Sachs-Wolfe plateau is a distinctive feature of Fig. These observations confirm the existence of a 
primordial spectrum of scalar (density) perturbations on all scales, otherwise the power spectrum would 
have started from zero all = 2. However, we see that the spectrum starts to rise around / = 20 towards 
the first acoustic peak, where the SW approximation breaks down and the above formulae are no longer 
vaUd. 

As mentioned above, the first peak in the photon distribution corresponds to overdensities that 
have undergone half an oscillation, that is, a compression, and appear at a scale associated with the size 
of the horizon at last scattering, about 1° projected in the sky today. Since photons scatter off baryons, 
they will also feel the acoustic wave and create a peak in the correlation function. The height of the peak 
is proportional to the amount of baryons: the larger the baryon content of the universe, the higher the 
peak. The position of the peak in the power spectrum depends on the geometrical size of the particle 
horizon at last scattering. Since photons travel along geodesies, the projected size of the causal horizon 
at decoupling depends on whether the universe is flat, open or closed. In a flat universe the geodesies 
are straight Unes and, by looking at the angular scale of the first acoustic peak, we would be measuring 



the actual size of the horizon at last scattering. In an open universe, the geodesies are inward-curved 
trajectories, and therefore the projected size on the sky appears smaller. In this case, the first acoustic 
peak should occur at higher multipoles or smaller angular scales. On the other hand, for a closed universe, 
the first peak occurs at smaller multipoles or larger angular scales. The dependence of the position of the 
first acoustic peak on the spatial curvature can be approximately given by [p^] 



^peak — 220 



-1/2 



(203) 



where Qq = Qm + = 1 — ^K- Present observations, specially the ones of the Mobile Anisotropy 
Telescope (MAT) in Cerro Tololo, Chile, which produced two data sets, TOC097 and TOC098 [ |ToT| ], 
and the recent balloon-borne experiment BOOMERANG [[7]]], suggest that the peak is between I = 180 
and 250 at 95% c.L, with an amplitude = 80 it 10 /xK, and therefore the universe is most probably 



flat, see Fig. 36, and Ref. [102]. In particular, these measuremts determine that 



0.85 < Sin < 1-25 



c.l. 



(204) 



That is, the universe is fiat, within 10% uncertainty, which is much better than we could ever do before. 
In the near future we will measure S7o to within 1%, with the new microwave anisotropy satellites. 
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Fig. 36: The left figure shows the power spectrum of the BOOMERANG experiment with 6 arcminute pixelization. The solid 
curve is a marginally closed model with {^Ib, ^m, Oa, n, h) — (0.05, 0.26, 0.75, 0.95, 0.7). The dotted curve is Standard 
CDM with (0.05, 0.95, 0.0, 1.0, 0.65). The dashed curves are open and closed models with fixed f^o = Qm + ~ 0.66 and 
1.55, respectively. The right figure shows the likelihood function of fio normalized to unity at the peak, after marginalizing 
over the Qm — direction. From Ref. |[7l|]. 

At the moment there is not enough information at small angular scales, or large multipole numbers, 
to determine the existence or not of the secondary acoustic peaks. These peaks should occur at harmonics 
of the first one, but are typically much lower because of Silk damping. Since the amplitude and position 
of the primary and secondary peaks are directly determined by the sound speed (and, hence, the equation 
of state) and by the geometry and expansion of the universe, they can be used as a powerful test of the 



density of baryons and dark matter, and other cosmological parameters, see Fig. |35. 

By looking at these patterns in the anisotropics of the microwave background, cosmologists can 



determine not only the cosmological parameters, see Fig. ^ but also the primordial spectrum of den- 
sity perturbations produced during inflation. It turns out that the observed temperature anisotropics are 



compatible with a scale-invariant spectrum, see Eq. ( |202| ), as predicted by inflation. This is remarkable, 
and gives very strong support to the idea that inflation may indeed be responsible for both the CMB 



anisotropies and the large-scale structure of the universe. Different models of inflation have different 
specific predictions for the fine details associated with the spectrum generated during inflation. It is these 
minute differences that will allow cosmologists to differentiate between alternative models of inflation 
and discard those that do not agree with observations. However, most importantly, perhaps, the pattern of 
anisotropies predicted by inflation is completely different from those predicted by alternative models of 
structure formation, like cosmic defects: strings, vortices, textures, etc. These are complicated networks 
of energy density concentrations left over from an early universe phase transition, analogous to the de- 
fects formed in the laboratory in certain kinds of liquid crystals when they go through a phase transition. 
The cosmological defects have spectral properties very different from those generated by inflation. That 
is why it is so important to launch more sensitive instruments, and with better angular resolution, to 
determine the properties of the CMB anisotropies. 




4.4.5 The new microwave anisotropy satellites, MAP and Planck 

The large amount of information encoded in the anisotropies of the microwave background is the rea- 
son why both NASA and the European Space Agency have decided to launch two independent satel- 
lites to measure the CMB temperature and polarization anisotropies to unprecendented accuracy. The 
Microwave Anisotropy Probe [72] will be launched by NASA at the end of 2000, and Planck [73] is 
expected in 2007. 

As we have emphasized before, the fact that these anisotropies have such a small amplitude al- 
low for an accurate calculation of the predicted anisotropies in linear perturbation theory. A particular 
cosmological model is characterized by a dozen or so parameters: the rate of expansion, the spatial cur- 
vature, the baryon content, the cold dark matter and neutrino contribution, the cosmological constant 
(vacuum energy), the reionization parameter (optical depth to the last scattering surface), and various 
primordial spectrum parameters like the amplitude and tilt of the adiabatic and isocurvature spectra, the 
amount of gravitational waves, non-Gaussian effects, etc. All these parameters can now be fed into a fast 
code called CMBFAST [103] that computes the predicted temperature and polarization anisotropies to 
1 % accuracy, and thus can be used to compare with observations. 

These two satellites will improve both the sensitivity, down to /iK, and the resolution, down to arc 
minutes, with respect to the previous COBE satellite, thanks to large numbers of microwave horns of var- 
ious sizes, positioned at specific angles, and also thanks to recent advances in detector technology, with 
high electron mobility transistor amplifiers (HEMTs) for frequencies below 100 GHz and bolometers for 
higher frequencies. The primary advantage of HEMTs is their ease of use and speed, with a typical sen- 
sitivity of 0.5 mKs^/^, while the advantage of bolometers is their tremendous sensitivity, better than 0.1 



see Ref. [|104|]. For instance, to appreciate the difference, compare the resolution in the temper- 
ature anisotropics that COBE and Planck would observe for the same simulated sky in Fig. This will 
allow cosmologists to extract information from around 3000 multipoles! Since most of the cosmological 
parameters have specific signatures in the height and position of the first few acoustic peaks, the higher 
the resolution, the more peaks one is expected to see, and thus the better the accuracy with which one 
will be able to measure those parameters, see Table 1 . As an example of the kind of data that these two 



satellites will be able to provide, see Fig. which compares the present observational status with that 
which will become available around 2008. 
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Fig. 38: The predicted angular power spectrum of temperature anisotropies, compared with the present data, see Fig. 
binned into 16 logarithmic intervals in multipole number between I — 2 and I — 1000. The right figure gives an estimate of the 
accuracy with which the power spectrum will be measured by Planck. It is only limited by cosmic variance on all the angular 



scales relevant to primary anisotropies. From Ref. [105] 



Although the satellite probes were designed for the accurate measurement of the CMB tempera- 
ture anisotropies, there are other experiments, like balloon-borne and ground interferometers, which will 
probably accomplish the same results with similar resolution (in the case of MAP), before the satellites 
start producing their own results [104]. Probably the most important objective of the future satellites will 
be the measurement of the CMB polarization anisotropies, yet to be discovered. These anisotropies are 
predicted by models of structure formation and are expected to arise at the level of micro Kelvin sensitivi- 
ties, where the new satellites are aiming at. The complementary information contained in the polarization 
anisotropies will provide much more stringent constraints on the cosmological parameters than from the 
temperature anisotropies alone. In particular, the curl-curl component of the polarization power spectra is 
nowadays the only means we have to determine the tensor (gravitational wave) contribution to the metric 
perturbations responsible for temperature anisotropies, see Fig. |3^. If such a component is found, one 
could constraint very precisely the model of inflation from its spectral properties, specially the tilt []100|]. 



4.5 From metric perturbations to large scale structure 

If inflation is responsible for the metric perturbations that gave rise to the temperature anisotropies ob- 
served in the microwave background, then the primordial spectrum of density inhomogeneities induced 
by the same metric perturbations should also be responsible for the present large scale structure [ ]106| ]. 
This simple connection allows for more stringent tests on the inflationary paradigm for the generation 
of metric perturbations, since it relates the large scales (of order the present horizon) with the smallest 
scales (on galaxy scales). This provides a very large lever arm for the determination of primordial spectra 
parameters like the tilt, the nature of the perturbations, whether adiabatic or isocurvature, the geometry 
of the universe, as well as its matter and energy content, whether CDM, HDM or mixed CHDM. 
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Fig. 39: Theoretical predictions for tiie four non-zero CMB temperature-polarization spectra as a function of multipole mo- 
ment. The dotted curves are the predictions for a COBE-normalized scalar perturbation from an inflationary model with no 
reionization and no gravitational waves for h = 0.65, VlW = 0.024, and A = 0. The solid curves are the corresponding 
predictions if the COBE anisotropy were entirely due to a stochastic gravitational wave background with a flat scale-invariant 
spectrum (with the same cosmological parameters). The panel for Cf^ contains no dotted curve because scalar perturbations 
produce no curl component of the polarization vector. From Ref. 1 100|. 



4. 5. 1 The galaxy power spectrum 

As metric perturbations enter the causal horizon during the radiation or matter era, they create density 
fluctuations via gravitational attraction of the potential wells. The density contrast 5 can be deduced from 



the Einstein equations in linear perturbation theory, see Eq. (134), 



5k = 

P 



- <Pfc = I I 



aH 



\ aH J 5 + 3a; 



(205) 



where we have assumed K = 0, and used Eq. (166). From this expression one can compute the power 



spectrum, at horizon crossing, of matter density perturbations induced by inflation, see Eq. (178), 



pik) = m' 



(206) 



with n given by the scalar tilt (180), n = 1 + 2?? — 6e. This spectrum reduces to a Harrison-Zel'dovich 
spectrum ( p9| ) in the slow-roll approximation: r?, e <C 1. 

Since perturbations evolve after entering the horizon, the power spectrum will not remain con- 
stant. For scales entering the horizon well after matter domination (k^^ » k^^ ~ 81 Mpc), the metric 



perturbation has not changed significantly, so that 7^^ (final) = 7?.^ (initial). Then Eq. (205) determines 
the final density contrast in terms of the initial one. On smaller scales, there is a linear transfer function 
T{k), which may be defined as 



7^fc(final) = r(/c) T^fc (initial) . 



(207) 



To calculate the transfer function one has to specify the initial condition with the relative abundance 
of photons, neutrinos, baryons and cold dark matter long before horizon crossing. The most natural 
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hot dark matter 




- 0.3 


5% 
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cosmological constant 




- 0.8 


8% 


0.5% 


spatial curvature 




0.2 - 1.5 


4% 


0.7% 


rate of expansion 
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0.4 - 0.8 
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age of the universe 
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spectral amplitude 
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20 - 30 


0.5% 


0.1% 


spectral tilt 
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0.5-1.5 


3% 


0.5% 


tensor-scalar ratio 
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- 1.0 


25% 


10% 


reionization 
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0.01 - 1.0 
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15% 



Table 1: The parameters of the standard cosmological model. The standard model of cosmology has around 12 different 
parameters, needed to describe the background space-time, the matter content and the spectrum of density perturbations. We 
include here the present range of the most relevant parameters, and the percentage error with which the microwave background 
probes MAP and Planck (without polarization) will be able to determine them in the near future. The rate of expansion is in 
units of Ho = 100 /i km/s/Mpc. 



condition is that the abundances of all particle species are uniform on comoving hypersurfaces (with 
constant total energy density). This is called the adiabatic condition, because entropy is conserved inde- 
pendently for each particle species X, i.e. 5px = Px^t, given a perturbation in time from a comoving 
hypersurface, so 

= , (208) 

PX +PX PY + PY 

where we have used the energy conservation equation for each species, px = —3H{px + Px)> valid to 
first order in perturbations. It follows that each species of radiation has a common density contrast 6r, 
and each species of matter has also a common density contrast 6m, with the relation 6m = j6r- 

Within the horizon, the density perturbation amplitude evolves according to the following equa- 



tion, see Ref. [|86p, 

H-% + [2 - 3(2^ - d)] H-'h-^il- Qcl + 8^ - 3a;^) 5, = - ' ^ 



(209) 



where oj = p/ p is the barotropic ratio, and (?s= pj p is the speed of sound of the fluid. 

Given the adiabatic condition, the transfer function is determined by the physical processes oc- 
curing between horizon entry and matter domination. If the radiation behaves like a perfect fluid, its 
density perturbation oscillates during this era, with decreasing amplitude. The matter density contrast 
living in this background does not grow appreciably before matter domination because it has negligible 



self-gravity. The transfer function is therefore given roughly by, see Eq. (52), 



1 , A; ^ fceq 



T{k) = 1 ^ " (210) 



The perfect fluid description of the radiation is far from being correct after horizon entry, because 
roughly half of the radiation consists of neutrinos whose perturbation rapidly disappears through free 
streeming. The photons are also not a perfect fluid because they diffuse significantly, for scales below 
the Silk scale, fc^^ ~ 1 Mpc. One might then consider the opposite assumption, that the radiation 



has zero perturbation after horizon entry. Then the matter density perturbation evolves according to 
Eq. ( |209| ), with 6 and p now referring to the matter alone, 



6k + 2H6k + {cl kl^ - 47rGp) 6^ = Q 



(211) 



which corresponds to the equation of a damped harmonic oscillator. The zero-frequency oscillator defines 
the Jeans wavenumber, kj = yJ'i-KGp/cj. For k <C kj, 6k grows exponentially on the dynamical 
timescale, Tdyn = Ima;~^ = (47rG/?)~^/^ = rgrav> which is the time scale for gravitational collapse. 
One can also define the Jeans length. 
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(212) 



which separates gravitationally stable from unstable modes. If we define the pressure response timescale 
as the size of the perturbation over the sound speed, Tprcs ~ A/c^, then, if Tprcs > 7"grav> gravitational 
collapse of a perturbation can occur before pressure forces can response to restore hydrostatic equilibrium 
(this occurs for A > Aj). On the other hand, if Tprcs < ''"grav. radiation pressure prevents gravitational 
collapse and there are damped acoustic oscillations (for A < A j). 

We will consider now the behaviour of modes within the horizon during the transition from the 
radiation (c^ = 1/3) to the matter era (c^ = 0). The growing and the decaying solutions of Eq. (ZH) are 



A[l + -y 
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(214) 



where A and B are constants, and y = a /acq. The growing mode solution ( |213| ) increases only by a 
factor of 2 between horizon entry and the epoch when matter starts to dominate, i.e. y = I. The transfer 



function is therefore again roughly given by Eq. ( |210| ). 

Since the radiation consists roughly half of neutrinos, which free streem, and half of photons, 
which either form a perfect fluid or just diffuse, neither the perfect fluid nor the free-streeming approx- 
imation looks very sensible. A more precise calculation is needed, including: neutrino free streeming 
around the epoch of horizon entry; the diffusion of photons around the same time, for scales below SiUc 
scale; the diffusion of baryons along with the photons, and the establishment after matter domination of a 
common matter density contrast, as the baryons fall into the potential wells of cold dark matter. All these 
effects apply separately, to first order in the perturbations, to each Fourier component, so that a linear 
transfer function is produced. There are several parametrizations in the literature, but the one which is 



more widely used is that of Ref. [107], 
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a = QA{VtMhy^ h'^ Mpc , 
6 = 3.0 (Om/i)"^ /i"^ Mpc, 
c = 1.7 {Quhy^ h-^ Mpc. 



(215) 

(216) 
(217) 
(218) 



We see that the behaviour estimated in Eq. (210) is roughly correct, although the break at A: = k^q is not 



at all sharp, see Fig. ^ The transfer function, which encodes the soltion to hnear equations, ceases to 
be valid when the density contrast becomes of order 1. After that, the highly nonlinear phenomenon of 



gravitational collapse takes place, see Fig. 40 
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Fig. 40: The CDM power spectrum P{k) as a function of wavenumber k, in logarithmic scale, normalized to the local abun- 
dance of galaxy clusters, for an Einstein-de Sitter universe with h = 0.5. The solid (dashed) curve shows the linear (non-linear) 
power spectrum. While the linear power spectrum falls off like k^^, the non-linear power-spectrum illustrates the increased 
power on small scales due to non-linear effects, at the expense of the large-scale structures. From Ref. ||47|]. 



4.5.2 The new redshift catalogs, 2dF and Sloan Digital Sky Survey 

Our view of the large-scale distribution of luminous objects in the universe has changed dramatically dur- 
ing the last 25 years [19|: from the simple pre-1975 picture of a distribution of field and cluster galaxies, 
to the discovery of the first single superstructures and voids, to the most recent results showing an almost 
regular web-Uke network of interconnected clusters, filaments and walls, separating huge nearly empty 
volumes. The increased efficiency of redshift surveys, made possible by the development of spectro- 
graphs and - specially in the last decade - by an enormous increase in multiplexing gain (i.e. the ability 
to collect spectra of several galaxies at once, thanks to fibre-optic spectrographs), has allowed us not 
only to do cartography of the nearby universe, but also to statistically characterize some of its properties, 
see Ref. [ ]109| ]. At the same time, advances in theoretical modeling of the development of structure, with 
large high-resolution gravitational simulations coupled to a deeper yet limited understanding of how to 
form galaxies within the dark matter halos, have provided a more realistic connection of the models to the 
observable quantities [ |11(]| ]. Despite the large uncertainties that still exist, this has transformed the study 
of cosmology and large-scale structure into a truly quantitative science, where theory and observations 
can progress side by side. 

For a review of the variety and details about the different existing redshift catalogs, see Ref. []T9|], 
and Fig. 41 . Here I will concentrate on two of the new catalogs, which are taking data at the moment and 
which will revolutionize the field, the 2-degree-Field (2dF) Catalog and the Sloan Digital Sky Survey 
(SDSS). The advantages of multi-object fibre spectroscopy have been pushed to the extreme with the 
construction of the 2dF spectrograph for the prime focus of the Anglo- Australian Telescope [111]. This 
instrument is able to accommodate 400 automatically positioned fibres over a 2 degree in diameter field. 
This implies a density of fibres on the sky of approximately 130 deg~^, and an optimal match to the 
galaxy counts for a magnitude bj ~ 19.5, similar to that of previous surveys like the ESP, with the 
difference that with such an area yield, the same number of redshifts as in the ESP survey can be collected 
in about 10 exposures, or slightly more than one night of telescope time with typical 1 hour exposures. 
This is the basis of the 2dF galaxy redshift survey. Its goal is to measure redshifts for more than 250,000 
galaxies with bj < 19.5. In addition, a faint redshift survey of 10,000 galaxies brighter than R = 21 
will be done over selected fields within the two main strips of the South and North Galactic Caps. The 
survey is steadily collecting redshifts, and there were about 93,000 galaxies measured by January 2000. 
See also Ref. [Ill], where the survey is continuously updated. 
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Fig. 41: Compilation of large-scale structure observations, showing the power spectrum P{k) as a function of wavenumber k. 
No corrections for bias, redshift distortions, or non-linear evolution have been made. Some of the redshift surveys have been 
rebinned to make the points nearly independent. The black box comes from measurements of erg from present-day number 
abundances of rich clusters, and the black point with error bars is from peculiar velocities. The height shows the 68% confidence 
interval. Q.m = 1 is assumed. The right panel shows a simulation of high-precision future CMB and LSS observations. MAP 
(red boxes) and Planck (blue boxes) are simulated assuming that CHDM is the correct model. Green error bars show the 
accuracy of the Sloan Digital Sky Survey and magenta error bars are for the 2 Degree Field Survey. No corrections are made 
for redshift distortions or non-linear evolution. The simulated data are indistinguishable from the underlying CHDM model for 
a wide range of k. From Ref. jlOq]. 



The most ambitious and comprehensive galaxy survey currently in progress is without any doubt 
the Sloan Digital Sky Survey [112]. The aim of the project is first of all to observe photometrically the 
whole Northern Galactic Cap, 30° away from the galactic plane (about 10^ deg^) in five bands, at limiting 
magnitudes from 20.8 to 23.3. The expectation is to detect around 50 million galaxies and around 10*^ 
star-like sources. This has already led to the discovery of several high-redshift {z > 4) quasars, including 
the highest-redshift quasar known, at 2 = 5.0, see Ref. [112]. Using two fibre spectrographs carrying 
320 fibres each, the spectroscopic part of the survey will then collect spectra from about 10^ galaxies 
with r' < 18 and 10^ AGNs with r' < 19. It will also select a sample of about 10^ red luminous galaxies 
with r' < 19.5, which will be observed spectroscopically, providing a nearly volume-limited sample 
of early-type galaxies with a median redshift of z ~ 0.5, that will be extremely valuable to study the 
evolution of clustering. The data expected to arise from these new catalogs is so outstanding that already 
cosmologists are making simulations and predicting what will be the scientific outcome of these surveys, 
together with the future CMB anisotropy probes, for the determination of the cosmological parameters 



of the standard model of cosmology, see Figs. ^ and ^ 



As often happens in particle physics, not always are observations from a single experiment suffi- 
cient to isolate and determine the precise value of the parameters of the standard model. We mentioned in 
the previous Section that some of the cosmological parameters created similar effects in the temperature 
anisotropics of the microwave background. We say that these parameters are degenerate with respect 
to the observations. However, often one finds combinations of various experiments/observations which 
break the degeneracy, for example by depending on a different combination of parameters. This is pre- 
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Fig. 42: Constraint regions in the Hm — Ho plane from various combinations of data sets. MAP data with polarization yields the 
ellipse from upper left to lower right; assuming the universe flat gives a smaller region (short-dashed line). SDSS (fcmax = O.lh 
Mpc^^) gives the vertical shaded region; combined with MAP gives the small filled ellipse. A projecton of future supernovae 
la results gives the solid vertical lines as bounds; combined with MAP gives the solid ellipse. A direct 10% measurement of 
Hq gives the long-dashed lines and ellipse. All regions are 68% confidence. The fiducial model is the f^M ~ 0.35 flat ACDM 



model. The right figure shows the same as before, but for constraints in the JIm — a plane. From Ref . [113] 



cisely the case with the cosmological parameters, as measured by a combination of large-scale structure 
observations, microwave background anisotropics, Supernovae la observations and Hubble Space Tele- 
scope measurements, a feature named somewhat idiosyncratically as "cosmic complementarity", see 



Ref. [113]. It is expected that in the near future we will be able to determine the parameters of the stan- 



dard cosmological model with great precision from a combination of several different experiments, as 



shown in Fig. 42 



5 CONCLUSION 

We have entered a new era in cosmology, were a host of high-precision measurements are already posing 
challenges to our understanding of the universe: the density of ordinary matter and the total amount 
of energy in the universe; the microwave background anisotropics on a fine-scale resolution; primordial 
deuterium abundance from quasar absorption lines; the acceleration parameter of the universe from high- 
redshift supernovae observations; the rate of expansion from gravitational lensing; large scale structure 
measurements of the distribution of galaxies and their evolution; and many more, which already put 



constraints on the parameter space of cosmological models, see Fig. gO[ However, these are only the 
forerunners of the precision era in cosmology that will dominate the new millennium, and will make 
cosmology a phenomenological science. 

It is important to bear in mind that all physical theories are approximations of reality that can fail 
if pushed too far. Physical science advances by incorporating earlier theories that are experimentally 
supported into larger, more encompassing frameworks. The standard Big Bang theory is supported by a 
wealth of evidence, nobody really doubts its validity anymore. However, in the last decade it has been 
incorporated into the larger picture of cosmological inflation, which has become the new standard cosmo- 
logical model. All cosmological issues are now formulated in the context of the inflationary paradigm. 
It is the best explanation we have at the moment for the increasing set of cosmological observations. 

In the next few years we will have an even larger set of high-quality observations that will test 
inflation and the cold dark matter paradigm of structure formation, and determine most of the 12 or 
more parameters of the standard cosmological model to a few percent accuracy (see table 1). It may 
seem that with such a large number of parameters one can fit almost anything. However, that is not 



the case when there is enough quantity and quality of data. An illustrative example is the standard 
model of particle physics, with around 21 parameters and a host of precise measurements from particle 
accelerators all over the world. This model is, nowadays, rigurously tested, and its parameters measured 
to a precision of better than 1% in some cases. It is clear that high-precision measurements will make the 
standard model of cosmology as robust as that of particle physics. In fact, it has been the technological 
advances of particle physics detectors that are mainly responsible for the burst of new data coming from 
cosmological observations. This is definitely a very healthy field, but there is still a lot to do. With 
the advent of better and larger precision experiments, cosmology is becoming a mature science, where 
speculation has given way to phenomenology. 

There are still many unanswered fundamental questions in this emerging picture of cosmology. 
For instance, we still do not know the nature of the inflaton field, is it some new fundamental scalar 
field in the electroweak symmetry breaking sector, or is it just some effective description of a more 
fundamental high energy interaction? Hopefully, in the near future, experiments in particle physics 
might give us a clue to its nature. Inflation had its original inspiration in the Higgs field, the scalar field 
supposed to be responsible for the masses of elementary particles (quarks and leptons) and the breaking 
of the electroweak symmetry. Such a field has not been found yet, and its discovery at the future particle 
colliders would help understand one of the truly fundamental problems in physics, the origin of masses. 
If the experiments discover something completely new and unexpected, it would automatically affect 
inflation at a fundamental level. 

One of the most difficult challenges that the new cosmology will have to face is understanding the 
origin of the cosmological constant, if indeed it is confirmed by independent sets of observations. Ever 
since Einstein introduced it as a way to counteract gravitational attraction, it has haunted cosmologists 
and particle physicists for decades. We still do not have a mechanism to explain its extraordinarily small 
value, 120 orders of magnitude below what is predicted by quantum physics. For several decades there 
has been the reasonable speculation that this fundamental problem may be related to the quantization of 
gravity. General relativity is a classical theory of space-time, and it has proved particularly difficult to 
construct a consistent quantum theory of gravity, since it involves fundamental issues like causality and 
the nature of space-time itself. 

The value of the cosmological constant predicted by quantum physics is related to our lack of 
understanding of gravity at the microscopic level. However, its effect is dominant at the very largest 
scales of clusters or superclusters of galaxies, on truly macroscopic scales. This hints at what is known 
in quantum theory as an anomaly, a quantum phenomenon relating both ultraviolet (microscopic) and 
infrared (macroscopic) divergences. We can speculate that perhaps general relativity is not the correct 
description of gravity on the very largest scales. In fact, it is only in the last few billion years that the 
observable universe has become large enough that these global effects could be noticeable. In its infancy, 
the universe was much smaller than it is now, and, presumably, general relativity gave a correct descrip- 
tion of its evolution, as confirmed by the successes of the standard Big Bang theory. As it expanded, 
larger and larger regions were encompassed, and, therefore, deviations from general relativity would 
slowly become important. It may well be that the recent determination of a cosmological constant from 
observations of supemovae at high redshifts is hinting at a fundamental misunderstanding of gravity on 
the very large scales. 

If this were indeed the case, we should expect that the new generation of precise cosmological 
observations will not only affect our cosmological model of the universe but also a more fundamental 
description of nature. 
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